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Abstract 

We derive a representation for the pion nucleon scattering ampli- 
tude that is vaUd to the fourth order of the chiral expansion. To 
obtain the correct analytic structure of the singularities in the low en- 
ergy region, we have performed the calculation in a relativistic frame- 
work (infrared regularization) . The result can be written in terms of 
functions of a single variable. We study the corresponding dispersion 
relations and discuss the problems encountered in the straightforward 
nonrelativistic expansion of the infrared singularities. As an applica- 
tion, we evaluate the corrections to the Goldberger-Treiman relation 
and to the low energy theorem that relates the value of the amplitude 
at the Cheng-Dashen point to the cr-term. While chiral symmetry does 
govern the behaviour of the amplitude in the vicinity of this point, the 
representation for the scattering amplitude is not accurate enough to 
use it for an extrapolation of the experimental data to the subthresh- 
old region. We propose to perform this extrapolation on the basis of a 
set of integral equations that interrelate the lowest partial waves and 
are analogous to the Roy equations for tttt scattering. 
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1 Introduction 



Here, perhaps an accuracy of 5 per cent is more appropriate. 

Henley and Thirring (||^], 1962) on the static model. 

The description of the pion nucleon interaction in terms of effective fields 
has quite a long history. The static model represents a forerunner of the 
effective theories used today. In this model, the kinetic energy of the nucleon 
is neglected: The nucleon is described as a fixed source that only carries spin 
and isospin degrees of freedom. For an excellent review of the model and its 
application to several processes of interest, we refer to the book of Henley 
and Thirring [|I|. 

The systematic formulation of the effective theory relies on an expansion 
of the effective Lagrangian in powers of derivatives and quark masses. Chiral 
symmetry implies that the leading term of this expansion is fully determined 
by the pion decay constant, F^, and by the nucleon matrix element of the 
axial charge, Qa- Disregarding vertices with three or more pions, the explicit 
expression for the leading term reads 

Ceff = -^i' 7^75 ^M^^ + 7''^[7r, df,'n]ip + ... 

The success of the static model derives from the fact that it properly accounts 
for the first term on the right hand side - in the nonrelativistic limit, where 
the momentum of the nucleons is neglected compared to the nucleon mass. 

The static model is only a model. In order for the effective theory to cor- 
rectly describe the properties of QCD at low energies, that framework must 
be extended, accounting for the second term in the above expression for the 
effective Lagrangian, for the vertices that contain three or more pion fields, 
for the contributions arising at higher orders of the derivative expansion, 
as well as for the chiral symmetry breaking terms generated by the quark 
masses m^. A first step in this direction was taken by Gasser, Sainio 
and Svarc 0, who formulated the effective theory in a manifestly Lorentz 
invariant manner. In that framework, chiral power counting poses a problem: 
The loop graphs in general start contributing at the same order as the cor- 
responding tree level diagrams, so that the loop contributions in general also 
renormalize the lower order couplings. A method that does preserve chiral 
power counting at the expense of manifest Lorentz invariance was proposed 
by Jenkins and Manohar 0, who used a nonrelativistic expansion for the nu- 
cleon kinematics. The resulting framework is called "Heavy Baryon Chiral 
Perturbation Theory" (HBCHPT). It represents an extension of the static 
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model that correctly accounts for nucleon recoil, order by order in the nonrel- 
ativistic expansion (for reviews of this approach, see for instance refs. [§, |^). 
In 10, 0], this formalism has been used to obtain an O(g^) representation 
for the pion nucleon scattering amplitude. In the meantime, the HBCHPT 
calculation has been extended to O^q"^) [§]. 

As pointed out in ref. 0, the nonrelativistic expansion of the infrared 
singularities generated by pion exchange is a subtle matter. The HBCHPT 
representations of the scattering amplitude or of the scalar nucleon form 
factor, for example, diverge in the vicinity of the point t = 4M^. The 
problem does not arise in the Lorentz invariant approach proposed earlier 
. It originates in the fact that for some of the graphs, the loop integration 
cannot be interchanged with the nonrelativistic expansion. 

The reformulation of the effective theory given in ref. exploits the fact 
that the infrared singular part of the one loop integrals can unambiguously 
be separated from the remainder. To any finite order of the nonrelativistic 
expansion, the regular part represents a polynomial in the momenta. More- 
over, the singular and regular pieces separately obey the Ward identities of 
chiral symmetry. This ensures that a suitable renormalization of the effective 
coupling constants removes the regular part altogether. The resulting rep- 
resentation for the various quantities of interest combines the virtues of the 
heavy baryon approach with those of the relativistic formulation of ref. 0: 
The perturbation series can be ordered with the standard chiral power count- 
ing and manifest Lorentz invariance is preserved at every stage of the calcu- 
lation. The method has recently been extended to the multi-nucleon sector 

Previously, this framework has been applied in calculations of the scalar, 
axial and electromagnetic form factors |^, 0, |12|. In the present paper, we 
use it to derive a representation for the pion nucleon scattering amplitude 
which is valid to the fourth order of the chiral expansion. After a compar- 
ison of infrared regularization with standard dimensional regularization, we 
discuss some technical issues in the evaluation of the Feynman diagrams. In 
sec. 1^ we give the chiral expansion of the masses and of the pion nucleon cou- 
pling constant g-„N- By comparing the quark mass dependence of the pion 
and the nucleon mass, we illustrate that the infrared singularities encoun- 
tered in the baryon sector are considerably stronger than in the Goldstone 
sector. 

In sees. p!0|-p!3| we discuss the constraints that chiral symmetry imposes on 
the scattering amplitude. Using the result for the axial coupling constant qa 
obtained by Kambor and Mojzis |TB|, we show that the Goldberger-Treiman 
relation is free of infrared singularities up to and including 0{q^). Our calcu- 



lation confirms a result obtained in fTl]: The difference between the E-Term 
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and the scalar form factor at t = 2M^ does not involve a chiral logarithm at 
order - the chiral expansion of the difference S — (t(2M^) starts with a 
term proportional to the square of the quark masses. 

In section |1^ we discuss the analytic structure of the amplitude. There 
are three categories of branch cuts in the one loop amplitude: cuts in the 
variables s and u due to vrA^-intermediate states and cuts arising from vrvr- 
and iVA^-intermediate states in the t-channel. The box graph is the only 
contribution to the amplitude which contains a simultaneous cut in two kine- 
matic variables. Since the NN cut only starts at t = 4m^, far outside the 
low energy region, its contribution to the amplitude is well represented by a 
polynomial. We make use of this fact to simplify the t-dependence of the box 
graph. Dropping higher order terms suppressed by t/4m^, the amplitude is 
given by the Born-term, a crossing symmetric polynomial of order and 9 
functions of a single variable, either s, t or u. These functions are given by 
dispersion integrals over the imaginary parts of the one loop graphs. 

The representation in terms of dispersion integrals contains terms of ar- 
bitrarily high order in the chiral expansion. The expansion of the integrals 
can be carried out explicitly. Truncating the series at O(g^), we obtain a 
simple, explicit representation of the scattering amplitude in terms of ele- 
mentary functions. The result can be compared directly with the representa- 
tion obtained in HBCHPT, where the chiral expansion of the loop integrals 
is performed ab initio. In contrast to that approach, our method allows us to 
examine the convergence of the chiral expansion of the loop integrals. The 



matter is discussed in detail in section |T8|, where we point out that the in- 
frared singularities contained in these integrals give rise to several problems. 
In particular, we show that the kinematic variables must be chosen carefully 
for the truncated expansion to represent a decent approximation. This leads 
to considerable complications in the analysis. Our dispersive representation 
of the scattering amplitude is comparatively simple and is perfectly adequate 
for numerical analysis. 

The representation of the scattering amplitude to 0{q^) neither includes 
the cuts due to the exchange of more than two stable particles nor the poles 
on unphysical sheets generated by resonances. It accounts for these effects 
only summarily, through their contribution to the effective coupling con- 
stants. Since the A-resonance lies close to the physical threshold it generates 
significant curvature in that region. In section |l^ we discuss the role of the 
higher order contributions both in the real and the imaginary part of the 
amplitude. 

Using the example of the S'-wave scattering lengths, we illustrate that the 
one loop representation of chiral perturbation theory does not cover a suffi- 
ciently large kinematic range to serve as a bridge between the experimentally 
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accessible region and the Cheng-Dashen point. The problem arises because 
the one loop amplitude fulfills the unitarity condition only up to higher or- 
ders in the chiral expansion and the perturbative series for the amplitude 
goes out of control immediately above threshold. In sec. ^ we discuss these 
unitarity violations on the basis of the optical theorem. 

We conclude that a reliable extrapolation from the physical region to the 
Cheng-Dashen point can only be obtained by means of dispersive methods 
and propose a set of integral equations for the lowest partial waves, similar 
to the Roy equations for Trvr-scattering. Based on partial wave relations 
IB] , p!7| , they allow for a determination of the S- and P-waves in the elastic 



region from their imaginary part at higher energies. The structure of the 
amplitude that underlies these equations matches the representation obtained 
in chiral perturbation theory, but in contrast to that representation, the 
lowest partial waves strictly obey the constraints imposed by unitarity. These 
equations allow us to extend the region where the one loop approximation 
of chiral perturbation theory provides a reliable description of the scattering 
amplitude. In the case of vrvr scattering, the method has been shown to 
yield a remarkably accurate representation of the amplitude throughout the 
elastic region |T^. We expect that the application of this method to nN 
scattering will lead to a reliable determination of the pion nucleon cr-term. 



2 Kinematics 




Figure 1: Kinematics of the elastic ttN scattering amplitude. P, q {P',q') 
denote the momentum of the incoming (outgoing) nucleons and pions, and 
a {a') stands for the isospin index of the incoming (outgoing) pion. 

We consider the scattering amplitude on the mass shell, 
P' = P" = ml, q' = q" = M'^. 



7 



In view of crossing symmetry, it is often convenient to express the energy 
dependence in terms of the variable 



s — u 

4mjv 

The standard decomposition involves the four Lorentz invariant amplitudes 

r± = v!{A^ + \{i' + i)B^]u. 

This decomposition is not suited to perform the low energy expansion of 
the amplitude, because the leading contributions from A and B cancel. We 
replace Ahy D = A + v B, so that the scattering amphtude takes the form 



u . 



Our evaluation of the chiral perturbation series to one loop allows us to 
calculate the amplitudes and B^ to 0{q^) and 0{q^), respectively. 



3 Effective Lagrangian 

The variables of the effective theory are the meson field U {x) e SU(2) and the 
Dirac spinor ^l){x) describing the degrees of freedom of the nucleon. The quark 
mass matrix = diag(m„, md) only enters the effective Lagrangian together 
with the constant B that specifies the magnitude of the quark condensate 
in the chiral limit. As usual, we denote the relevant product by % = 2Bmg. 
The constant B also determines the leading term in the chiral expansion of 
the square of the pion mass, which we denote by M^, 

= (m, + md)B . 

We disregard isospin breaking effects and set m„ = m^. 

The effective Lagrangian relevant for the sector with baryon number equal 
to one contains two parts. 

The first one only involves the field U (x) and an even number of derivatives 
thereof. Using the abbreviations 
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this part of the Lagrangian takes the form [| 

= lF2(«X + X+) + |^4(«X)(x+) + i^(^3 + /4)(x+)' + 0(g^). (3.1) 

The expression agrees with the one used in ref. 0, but differs from the La- 
grangian in ref. pO[ by a term proportional to the equation of motion. The 
difference is irrelevant, because it amounts to a change of the meson field 
variables. For the physical quantities to remain the same, this change of 
variables, however, also needs to be performed in the term Cm - the signifi- 
cance of the effective coupling constants occurring in that term does depend 
on the specific form used for £^ pT| . 

The term is bilinear in ip{x), ijj{x) and the derivative expansion con- 
tains odd as well as even powers of momentum: 

The leading contribution is fully determined by the nucleon mass m^v, the 
pion decay constant F^^ and the matrix element of the axial charge qa- We 
denote the values of these quantities in the chiral limit by m, F and g, 
respectively. With = d^ + T^, the explicit expression then takes the form 

The Lagrangian of order contains four independent coupling constants^ 
= ci ^ - ^ (u^u,) ii' D^D-"^ + h.c.) (3.2) 

The terms and C^^^ enter the representation of the scattering amplitude 
to 0{q^) only at tree level. The Lagrangian of order 0{q^) is taken from 

^In the notation of ref. |2^, the matrix U stands for U = U'^+i f-U. In the isospin Hmit, 
we have x+ = 2M^U°, x-=-^^M^t-U, so that (x-) = 0, (xl) = - Hxh)- 

^We use the conventions of ref. ||2^. In this notation, the couphng constants of ref. [|| 
are given by: m cf^^ = F^ci , m c§^^ = — F^C4, m c^^^ = — 2F^C3, m (cf^'^ —2m c^^^) = 
— F^C2 (to order q^, the terms cf^^ and c^^^ enter the observables only in this combination), 
while those of ref. pl[ read: 16 ai — Smc^ + g^, 8a2 — 4,mc2 — g^, 03 = mci, 4a5 = 
4mc4 + 1 — g^. In the numerical analysis, we work with F^^ = 92.4MeV, gA — 1.267, 
mjv ~ mp, = M^+. 
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ref. [§: 

+ y|^(K, [Du, ux]]iD^D''D' + sym.) + h.c.) 
+ ^([X-,-.]^'^ + h.c.) 

At tree level only five combinations of the coupling constants occurring here 
enter the scattering amplitude. We do not give an explicit expression for 
C^^\ but characterize this part of the effective Lagrangian through its con- 
tributions to the scattering amplitude (sec. ^). 

4 Tree graphs 

The tree graph contributions to the amplitude have the following structure 




The vertex contributions (s) involve the axial coupling g and a quark mass 
correction to it from the Lagrangian C^^^ . The quark mass correction can 
be taken into account by replacing the coupling constant g in C^^^ with 
92 = 9 + 2M^(2(ii6 — dis). The contributions of type @ represent mass 
insertions generated by and at tree level. They are taken into 
account by replacing the bare mass m in with = m — AciM^ + ei 
(we adopt the notation of and denote the relevant coupling constant in 
by ei). 

When evaluating the relevant diagrams, we must distinguish between the 
bare and physical masses of the nucleon. As we are working on the physical 
mass shell, = m^, we need to use the corresponding Dirac equation 
for the nucleon spinors, f u{P) = mf^u{P). The result for the sum of all 
contributions of type@ and ©reads (here, the amplitude A is more convenient 
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than D, because the representation then only involves functions of a single 
variable): 

{rrii + m^r) f s - ml 



A^=A{s)±A{u) A{s) 



4F2 



= Bis) T Bin) Bis) = + + ^ 



4 F2 (ml - s) 

The tree graphs of type @ yield polynomials in z/, t and M'^. We only list 
the nonzero contributions. Those from are given by 

V 1 

D- = —r , B 



2p2 ' 2F2 ' 

while generates the terms 

4ciM2 ^ C2(16m>2_t2) ^ c3(2M2-t) ^__2c4m;, 



Note that the coupling constant C2 gives rise to a contribution proportional 
to which is of 0{q^). Quite generally, vertices involving derivatives of the 
nucleon field may generate contributions beyond the order indicated by chiral 
power counting. 

In the graphs from and dt\ we may replace the bare constants 
by the physical ones, because the distinction is beyond the accuracy of our 
calculation. The term C)^ only shows up in those amplitudes that are odd 
under crossing: 

D- = {2(c?i + d2 + 2^5) Ml - (rfi + d2)t + 2d^ u^} + Oiq') 

n 

B+ = ^^{d,,-d,,)+0{q') 

n 

Finally, generates a polynomial that is even under crossing. We identify 
the coupling constants 63, ... , en with the coefficients of this polynomial: 

D+ = l^{e,M^ + e^My + e,Mlt + e,u' + e7uH + est^}+0{q^) 

= ^ {e,Ml + eio ly' + en t} + 0{q') (4.1) 

A contribution analogous to e^M^ + e^M^ t also occurs in the representation 
of the scalar form factor to O(g^) (in P] the corresponding contribution to 
ait) is denoted by 2ei + e2Mlt). The coupling constants are different, 
because the terms and {x+){Uf_iU'^) ijjil' oi C^^^ do contribute to the 

scattering amplitude, but do not show up in the scalar form factor. 
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(b) 
-t r 



(c) 




(n) (o) (p) 




W (s) (t) 




Figure 2: One loop topologies for TrA^-scattering. Crossed diagrams and 
external leg corrections are not shown. We do not display topologies involving 
closed fermion loops since they vanish in infrared regularization. 
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5 One loop graphs in infrared regularization 



To obtain a representation of the scattering amplitude to order q^, it suffices 
to evaluate the one loop graphs of . As we wish to evaluate the scattering 
amplitude to order g^, we in addition need to consider loops containing one 
vertex from C)^ ■ The various topologies are shown in figure In this figure, 
we do not differentiate between the vertices of £^'* and those of . Graph 
(a) for instance describes two different contributions (see figure]^). Since the 




0{q') 0{q') 

Figure 3: By exclusively inserting vertices from (drawn as full circles) 
into topology (a), we obtain a graph of O(g^). The same topology also 

(2) 

occurs if one of the vertices is replaced by one of those in , which involve 
the coupling constants Ci, C2, C3 or C4 (denoted by a circled cross). The 
resulting graph then starts contributing at O(g^). Since exclusively 
contains vertices with an even number of pions, only the rightmost vertex 
can be replaced. 

second order Lagrangian C)^ only involves couplings with an even number 
of pions, the topologies (c), (d), (e), (i), (1), (p), (r), (t) and (u) do not give 
rise to 0{q'^) diagrams. 

The calculation of the loop diagrams generated by was carried out 
quite some time ago by Gasser, Sainio and Svarc 0, using a relativistic La- 
grangian and dimensional regularization. As mentioned in the introduction, 
this approach does not preserve the counting rules. In a recent paper , we 
have have set up a variant of dimensional regularization that avoids these dif- 
ficulties and which we call "infrared regularization" . We now briefly discuss 
this method. 

5.1 Infrared regularization 

To explain the essence of the method, we consider the simplest example, the 
self energy graph shown in figure The corresponding scalar loop integral 
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k 



■> 



P 



P-k 



P 



Figure 4: Self energy 



has the form 



H 



d'^k 1 1 




d'^k 1 1 
(27rp a h 



(5.1) 



{2'kY M2 - F m2 - (P - k) 



We need to analyze the integral for external momenta in the vicinity of the 
mass shell: P — mv + r, where v is a timelike unit vector and r is a quantity 
of order q. In the limit M — > 0, the integral develops an infrared singularity, 
generated by small values of the variable of integration, k = 0{q). In that 
region, the first factor in the denominator is of 0{q'^), while the second is of 
order 0{q). Accordingly the chiral expansion of the integral contains terms of 
order q'^~^. The order of the infrared singular part follows from the counting 
rules at tree level, because it is generated by the region of integration, where 
the momenta flowing through the propagators arc of the same order as in the 
tree level graphs. The remainder of the integration region does not contain 
infrared singularities and may thus be expanded in an ordinary Taylor series. 
An evaluation of the integral at = s+ = (m + M)^ nicely shows the two 
parts: 



The infrared singular part / is proportional to M'^~^, while the remainder 
R is proportional to m'^~^ and does therefore not contain a singularity at 
M = 0, irrespective of the value of d. Since the regular contribution stems 
from a region where the variable of integration k is of order m, it violates 
the tree-level counting rules. 

The explicit expression for the loop integral H[P^) in d dimensions is 
lengthy, but the splitting into the infrared singular part / and regular part R 
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is easily obtained in the Schwinger-Feynman representation of the integraL 



Schwinger-Feynman parametrizations of the infrared singular and regular 
parts of a general one loop integral are given in 0. There, the following 
statements are proven: 

• To any given order, the chiral expansion of the infrared regular part of 
a one loop integral is a polynomial in the quark masses and the external 
momenta. 

• The regular part of the one loop amplitudes is chirally symmetric. 

Taken together, the two statements imply that, in a one loop calculation, 
one may simply drop the regular parts of the loop integrals because the 
contributions they generate can be absorbed in the coupling constants of the 
effective Lagrangian. This procedure is called "infrared regularization" : The 
loop graph contributions are identified with the infrared singular parts of the 
corresponding dimensionally regularized integrals. 

5.2 Comparison with dimensional regularization 

The calculation of loop graphs in infrared regularization is completely anal- 
ogous to the evaluation in dimensional regularization. The expressions for 
the graphs in terms of loop integrals look the same. In dimensional regu- 
larization, the contribution of the graph (f) in figure ^ to the isospin even 
amplitude, for instance, is given by 



The terms A^r and A^r stand for the pion and nucleon propagators at the 
origin, respectively: 
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The function H^^\s) may be expressed in terms of the scalar self energy 
integral H{s), which we considered in the preceding subsection: 

H^'\s) = 1 {(s-m' + M') H{s) + A. - A^} . 

i-i s 

To obtain the corresponding result in infrared regularization, it suffices 
to replace all loop integrals by their infrared singular parts. Integrals like A^ 
with pion propagators only do not have an infrared regular part, because they 
do not involve the heavy scale m: In the meson sector, infrared regularization 
coincides with ordinary dimensional regularization. Nucleon loop integrals 
like Ajv, on the other hand, vanish in infrared regularization, because they are 
infrared regular. Integrals which involve both kinds of propagators do posses 
a regular and an infrared singular part. These integrals are replaced by their 
infrared singular part. The result for graph (f) in infrared regularization 
is thus obtained with the substitution H(s) — > /(s), H^^\s) I^^\s) and 
A^ ^0. 

The result for the various graphs is given in appendix 0. The definitions 
of the pertinent loop integrals can be found in appendix O. We have checked 
that our result for the loop graphs of C]^ agrees with the one of Gasser, Sainio 
and Svarc 0, if we evaluate our loop integrals in dimensional regularization 
and expand around (1 = 4. Note that we do not display the contributions 
proportional to loop integrals that exclusively contain nucleon propagators, 
because these vanish in infrared regularization. Also, the basis used in the 
decomposition of the vector and tensor integrals differs from the one used by 
Gasser, Sainio and Svarc - in our basis, the chiral power counting is more 
transparent. 

6 Simplification of the 0{q^) diagrams 

The chiral expansion of the loop graphs in general contains terms of arbitrar- 
ily high order. Since we need the amplitude only to 0(g'^), we can simplify 
the representation by neglecting terms of 0{q^) or higher. The example of 
the triangle graph (graph h in fig. |^) shows, however, that the chiral expan- 
sion of infrared singularities is a subtle matter ^ , which we will discuss in 
detail later on. The problem arises from the denominators that occur in the 
various loop integrals - the expansion thereof gives rise to an infinite series 
of terms, which in some cases needs to be summed up to all orders to arrive 
at a decent representation of the integral. Concerning, the numerators, this 
problem does not arise - the corresponding expansion does not give rise to 
an infinite series of terms and hence preserves the analytic structure of the 
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integral, term by term. We exploit this fact to simplify the calculation of 
those graphs that contain vertices from C)^ , which start contributing only 
at O(g^) (in the loop graphs generated by C^k\ we retain all contributions). 

As an illustration, we consider the part of the 0{q'^) graph (a) that is 
proportional to C2 

_ C29I 1 I f d^k 1 



4F4m2 m^-s i j i {2tiY {M^ - k^){m^ - {P + q - k)^) 
X {k-P'q'-P' + k-{P + q-k) q'-{P + q-k)) 

X u' (m + f+^-^)75^(m + /'+^)75^M 

The numerator of the integrand involves terms with up to five powers of the 
loop momentum. It can explicitly be expressed in terms of the functions 
/(s), that arise in the tensorial decomposition of the generic 

loop integral with one meson and one nucleon propagator. The resulting 
expression is rather lengthy, however. As just discussed, we only retain the 
leading term in the expansion of the numerator: 

_ C29I 1 1 M'^ 



x2k-P'q ■ P'u2m-f5lt{m + f ) 75^^ + 0{q^) 
02 9 A 1 /■ d'^k 1 



« J I 



2F^m i Ji {2t:Y (M^ - k"^) (m'^ - (p + q - kf) 

X k- P'u^m- f)iu + 0(q^) , 

where we have used P' = P + 0{q). The relation 2k- P' = rn? — {P + q — k)"^ + 
0{q^) then allows us to remove the nucleon propagator, so that we arrive at 

The remaining integral vanishes because the integrand is antisymmetric in 
k, so that the entire contribution is of 0{q^) and can be dropped. 

7 Wave function renormalization 

According to the rules of perturbation theory, the scattering amplitude may 
be evaluated from the connected four-point-function 

{0\TN{x')Pa>{y')Pa{y)N{x)\0). 
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where N{x) and Pa{x) are any operators that interpolate between the relevant 
incoming and outgoing states, so that the matrix elements 

{0\N{x)\N) = yz^ue-^^^ , {0\Pa{x)\7rt) = ^.ba^e^'^^ , 

are different from zero (we suppress the spin and isospin quantum numbers 
of the nucleon). The scattering amplitude T^'a is determined by the residue 
of the poles occurring in the Fourier transform of this correlation function: 



t d^x'dWd'y{0\TNix') Paiy') P„(y) Nix) |0)^e 



i(P'x'+q'y'-Px-qy) 



U Ula'a + ■ ■ ■ 



Ml - g'2 M2 - g2 ^2 _ p/2 ^2 _ p2 

The result for Ta'a does not depend on the choice of the interpolating fields. 

In our context, it is convenient to identify the interpolating fields with the 
variables of the effective theory, N{x) = ipi^x), Pa{x) = Ha{x). We emphasize 
that these are auxiliary quantities that represent the variables of integration 
in the path integral - the effective fields do not have counterparts in QCD 
and do thus not represent objects of physical significance. The renormal- 
ization of the coupling constants occurring in our effective Lagrangian only 
ensures that the scattering amplitude remains finite when the regularization 
is removed, d —>■ 4. A multiplicative renormalization of the effective fields, 
ipix) = ipixY"'', TCaix) = ^fZ^TiJyxY'''^ ■, docs uot sufficc to obtain fi- 
nite correlation functions for these. Note also that the correlation functions 
depend on the parametrization used for the matrix U[x). 

It is easy to see why the correlation functions of the effective fields do not 
represent meaningful quantities. The corresponding generating functional is 
obtained by adding source terms to the effective Lagrangian, 

This operation, however, ruins the symmetries of the effective theory, because 
the representation of the chiral group on Tia{x) and il){x) involves the pion 
field in a nonlinear manner. To arrive at finite correlation functions, we 
would need to add extra counter terms that are not invariant under chiral 
rotations. 

A more physical choice for the interpolating fields is discussed in ap- 
pendix |B[ This is of interest, for instance, in connection with the matrix 
elements relevant for baryon decay. The corresponding correlation functions 
do have a physical interpretation within QCD. Since the relevant operators 
carry anomalous dimension, they depend on the running scale of that theory. 
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but are otherwise free from ambiguities. In our context, the choice of the 
interpolating fields is irrelevant - the result obtained for the S-matrix by 
using the fields discussed in appendix ^ would be the same. 

The perturbative calculation of the scattering matrix can be simplified 
in the familiar manner: It suffices to consider the amputated four-point- 
function, obtained by (i) discarding all graphs that contain insertions in 
the external lines and (ii) replacing the free propagators that describe the 
external lines in the remaining graphs by plane wave factors. Denoting the 
amputated four-point-function by Ta'a, the scattering amplitude is given by 

Ta'a = Zj^ Z^ U Ta'a U ■ 

We have performed our calculations in the so-called sigma parametrization 



TTlXl^ . Tf(x) ■ T 



U{x) = \ 



F2 F ' 

The wave function renormalization constants are then given by 

2 

Zn = l-^{h + K) + 0{M') 



A. 



2F2 [ 487r2 327rm 

M"^-^ r 1 1 



(^ln47r + r(l) + l)| 



(47r)2 [d-A 2 

The divergences oc (d — 4)~^ are contained in the quantity A^^, which involves 
the pion mass instead of an arbitrary scale /i - we are in effect identifying 
the renormalization scale with the pion mass. This is convenient, because it 
simplifies the formulae: The divergences are always accompanied by a chiral 
logarithm. Expressed in terms of the usual pole term A, we have 

The effective coupling constants of C^^\ and C^^^ pick up renormaliza- 
tion. The relevant formulae are listed in appendix |^. We have checked 
that if the scattering amplitude is expressed in the renormalized couplings 
il{fi),il{^),d'-{fi),el{fi) introduced there, it indeed remains finite at — > 4. 

8 Chiral expansion of niN-, Qtzn and 

The graphs (e) and (s) in fig. |^ contain a double pole at s = m\ and their 
crossed versions exhibit the analogous singularity at m = m\ (we recall that 
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the mass shifts generated by the tree graphs of and are included 
in the free Lagrangian, so that the nucleons propagating in the loops are 
equipped with the mass m4 = m — 4ciM2 + ei M^). The double pole dis- 
appears as it should if the Born term, which contains the pole due to one- 
nucleon exchange, is expressed in terms of the physical mass of the nucleon. 
In the context of the low energy expansion, the standard, pseudoscalar form 
of the Born term is not convenient, because it corresponds to a scattering 
amplitude of order We instead use the pseudovector form, defined by 

^^.(^,t) = ^^^, i?>,t) = ^V^, (8.1) 

D7,M. t) = uB-(^, t) , B-(u, t) - 



2ml 



s-u t- 2M: 



2 



^ = -; 1 



4 m„ 4 m. 



The Born term exclusively involves the physical values of m^r, g-KN and M^^. 
The explicit expression for rrii^ reads ||^ 

, , 3 o^M^ 3 (o^ - 8 ci m + C2 m + 4 cg m) A^M^ 
= "-4c,Af 

^<2»'-^-'")*^\e,M' + 0(M=). (8.2) 



1287r2 m 

while the one for the pion-nucleon coupling constant is given by 

mg r 4ciM2 2 (2 die - dis) ^g'^KM'^ 

fi-vriV = ^ i 1 



F \ F"^ m g F^ 

g'M' ^ il2 + 3g'-lQcsm + 32c,m)M' ^ ^^^^^,^\ ^ ^^^^^ 



167r2F2 967rmF2 

To order , the result agrees with the one given in ref . - only the terms of 
order are new. The quantity g represents the value of the axial charge in 



the chiral limit. Kambor and Mojzis [0] and Schweizer |TT| have calculated 
the physical axial charge gA to 0{q^), with the result 

, 4rfi6M2 2{2g^ + l)KM^ g^ 

9a = g<l + 



g F2 167r2F2 
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These expressions are analogous to the well-known representations that 
describe the dependence of M^, and F^, on the quark masses [pO|] , 



Ml = M^h + — {2h + K) + 0{M') 

F^ = f|i + — (/4-2A^) + 0(M^: 

In all cases, the expansion contains contributions that are not analytic in the 
quark masses, generated by the infrared singularities. In the meson sector, 
only the chiral logarithm contained in A^r occurs, but in the nucleon sector, 
the expansion in addition involves contributions with oc (m„ + m^) 2 . 



9 Strength of infrared singularities 

The fact that the low energy expansion in the baryon sector contains even 
as well as odd powers of M implies that the infrared singularities occurring 
there are stronger than those in the meson sector. We illustrate this with 
the (T-term matrix element, that is, with the response of the nucleon mass to 
a change in the quark masses: 

(y = mu- l-rud- . (9.1) 

Oniu OTTld 

More specifically, we wish to discuss the dependence of this quantity on the 
quark masses, which is quite remarkable. In this context, the numerical value 
of a is not of crucial importance. For definiteness, we use a = 45 MeV ||23[ . 
The chiral expansion of the a-term is readily obtained from the formula ( ^.2|) , 
with the result 

a = kiM^ + - k2M^ + k^M^ <^ 2 In — + 1 ^ + 2 k^M^ + O(M^) , (9.2) 
2 y ] 

As discussed above, the term proportional to arises from an infrared 
singularity in the self energy of the pion cloud. It lowers the magnitude 
of the cr-term by 3/2 x 15 MeV ~ 23 MeV. The coefficient k^ can also be 
expressed in terms of measurable quantities Numerically the contribution 
from this term amounts to —7 MeV, thus amplifying the effect seen at O(M^). 
Chiral symmetry does not determine all of the effective coupling constants 
entering the regular contribution k^M'^, which is of the same type as the 
correction = kcoM^ to the low energy theorem (|11.3|) . As discussed 
above, corrections of this type are expected to be very small - we simply 
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T 

Figure 5: cr-term (in MeV) as a function of the quark masses. The hori- 
zontal axis gives the values of = (m„ + md)B. The arrow corresponds to 
the physical value of the quark masses. For definiteness, it is assumed that 
the physical value of a is 45 MeV (dash-dotted line). The dashed hne de- 
picts the linear dependence that results if the infrared singular contributions 
proportional to and to In M^/m^ are dropped. 



drop the term k^M^. The value of ki is then fixed by the input cr = 45 MeV 
for the total, so that we can now discuss the manner in which a changes 
when the quark masses are varied. 

At leading order, a is given by kiM^. In figure 1 this contribution is 
shown as a dashed straight line. The full curve includes the contributions 
generated by the infrared singularities, k2M^ and A;3M^{2 In M^/m^ + 1}. 
The figure shows that the expansion of the cr-term in powers of the quark 
masses contains large contributions from infrared singularities. These must 
show up in evaluations of the cr-term on a lattice: The ratio a/{mu + ma) 
must change significantly if the quark masses are varied from the chiral limit 
to their physical values. Note that in this discussion, the mass of the strange 
quark mass is kept fixed at its physical value - the curvature seen in the 
figure exclusively arises from the perturbations generated by the two lightest 
quark masses. 

It is instructive to compare this result with the dependence of on the 
quark masses. In that case, the expansion only contains even powers of M: 

^ + 32^ AT + ^^^'^ ■ ^^-^^ 

The quantity A3 stands for the renormalization group invariant scale of the 
effective couphng constant ^3. The SU(3) estimate for this coupling constant 
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Figure 6: Square of the pion mass as a function of the quark masses. The 
horizontal axis gives the values of the leading order term = {mu + md)B. 
The dash-dotted line indicates the physical value of M^. 



given in ref. |2J] reads 1^ = — InM^/Ag = 2.9 ± 2.4. The error bar is so 
large that the estimate barely determines the order of magnitude of the scale 
A3. Figure 2 shows, however, that this uncertainty does not significantly 
affect the dependence of on the quark masses, because the logarithmic 
contribution is tiny: the range of I3 just quoted corresponds to the shaded 
region shown in the figure. 

The logarithmic term occurring in the chiral expansion of Ml gets en- 
hanced by about a factor of two if we consider the pion a-term, 

(9M2 dM'^ 

cr^ = (7r| ruu uu + dd \v:) = ruuirr-^ + rrid^r—^ . (9.4) 

omu onid 

We do not show the corresponding curve, because it is also nearly a straight 
line. 

The main point here is that the infrared singularities encountered in the 
self energy of the nucleon are much stronger than those occurring in the self 
energy of the pion. 



10 Goldberger-Treiman relation 

The above expressions for the quantities t^^jv, gj,, rrit^ and allow us to 
analyze the correction that occurs in the Goldberger-Treiman relation up to 
and including terms of order M^. We write the relation in the form 

= ^^{1 + Aar} ■ (10.1) 
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If the quark masses niu and are turned off, tfie strength of the vrA^ inter- 
action is fully determined by and F^: Aqj, = 0. Inserting the formulae 
of the preceding section in the ratio g/rrij^/FT, and expanding the result to 
order M^, we find that the coefficients of the infrared singular terms pro- 
portional to In M^/m^ and to are identical with those in (7^^^. Up 
to and including third order, the correction Aqt is therefore free of infrared 
singularities and takes the simple form 

^ - ^^^^^'+0(M^). (10.2) 



9 

As far as the logarithmic terms of O(M^) are concerned, this result was 
established a long time ago 0. What the present calculation adds is that 
the infrared singularities also cancel at O(M^), in agreement with ref. 

The result shows that in the case of the Goldberger-Treiman relation, 
the breaking of chiral symmetry generated by the quark masses does not get 
enhanced by small energy denominators. Chiral symmetry does not deter- 
mine the magnitude of the coupling constant dis- A crude estimate may 
be obtained from the assumption that the scale of the symmetry breaking 
is the same as in the case of Fk/Ft^, where it is set by the massive scalar 
states, Ms ~ 1 GeV. This leads to A^t ^ Ml/Mj ~ 0.02. A more de- 
tailed analysis based on models and on the SU(3) breaking effects seen in 



the meson-baryon coupling constants may be found in ref. [^. Recently 
the Goldberger-Treiman discrepancy has also been determined from QCD 
sum rules ||2^. Both evaluations confirm the expectation that A^t niust be 
very small - a discrepancy of order 4% or more would be very difficult to 
understand. 



Since the days when the Goldberger-Treiman relation was discovered P7| , 
the value of Qa has increased considerably. Also, decreased a little, on 
account of radiative corrections. The main source of uncertainty is q-j^m- 

The 



comprehensive analysis of tiN scattering published by Hohler in 1983 |2^ led 
to p = gl^Mll ilQ-wml) = 0.079. With = 1.267 and F^ = 92.4 MeV, this 
value yields Aqt = 0.041. The data accumulated since then indicate that 
is somewhat smaller, numbers in the range from 0.076 to 0.077 looking more 
likely. This range corresponds to 0.021 < Aar < 0.028. 

We conclude that, within the current experimental uncertainties to be 
attached to the pion-nucleon coupling constant, the Goldberger-Treiman re- 
lation does hold to the expected accuracy. Note that at the level of 1 or 2 %, 
isospin breaking cannot be ignored. In particular, radiative corrections need 
to be analyzed carefully. Also, the coupling constant relevant for the neutral 
pion picks up a significant contribution from vr" — 77 interference. 
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11 Low energy theorems for and D 

The representation of the scattering amphtude to order also allows us 
to analyze the corrections occurring in the well-known low energy theorem 
PU| , Pm for the value of the scattering amplitude D~^{h',t) at the Cheng- 
Dashen point, u = 0, t = 2M^. The theorem relates the quantity^ 

E = F25+(0,2M2) (11.1) 

to the scalar form factor of the nucleon, 

{N'\muUU + mddd\N) = a{t)U'u . (11.2) 

The relation may be written in the form 

S = (t(2M2) + Aez, . (11.3) 

The theorem states that the term Aqd vanishes up to and including contri- 
butions of order M^. The explicit expression obtained for S when evaluating 
the scattering amplitude to order again contains infrared singularities pro- 
portional to and In M"^ /m? . Precisely the same singularities, however, 
also show up in the scalar form factor at t = 2M^, so that the result for Ac^, 
is free of such singular it ies0: 

Acz, = f-^T^ - 2 ei - 2 62 + 63 + 2 65 + 4 eg) + 0{M^) . (11.4) 

Equally well, we could replace the bare coupling constants by the renor- 
malized ones: In the combination occurring here, the poles at c? = 4 con- 
tained in 6i, . . . ,68 cancel. Crude estimates like those used in the case of the 
Goldberger-Treiman relation indicate that the term A^o must be very small, 
of order 1 MeV. 

Unfortunately, the experimental situation concerning the magnitude of 
at the Cheng-Dashen point leaves much to be desired (for a recent review, 
see ref. [Q). The low energy theorem makes it evident that we are dealing 
with a small quantity here - the object vanishes in the chiral limit. The 
inconsistencies among the various data sets available at low energies need to 
be clarified to arrive at a reliable value for (7^^. Only then will it become 
possible to accurately measure small quantities such as the cr-term. 
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At the Cheng-Dashen-point, the amphtude is singular, on account of the Born term 



The bar indicates that this term is removed in the pseudovector form specified in eq. (B.l) 
= ~ = - B^ 

"'The cancellation of the terms of order M'^ was pointed out in ref. ^ and the 
absence of logarithmic contributions of order was shown in ref. . 
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There is a low energy theorem also for the isospin odd amplitude: In the 
chiral limit, the quantity 



C = 2F: 



D-{u,t) 



= 0, t = 2Ml 



:il.5) 



tends to 1. The representation of the scattering amplitude to O(g^) yields 
the corrections up to and including O(M^): 



C 
ki 
k2 



1 - 



M 



— — - M2 In— + kiM^ + k2 + 0{M^) 



;ii.6) 



;i4-3 7r) + (22 + 3 7r)^; 



2887r2F2 



64mci + 2 gl{A + gl) + ^/2 ln(l + V2) g 



Note that the scale dependence of the renormalized coupling constant d^fi) 
cancels against the one of the chiral logarithm. 

In contrast to the quantity Aqj^, this expression does contain infrared 
singularities at the order considered here: A chiral logarithm as well as terms 
of order M^. We will discuss the size of these effects in section O- first, we 



need to establish contact with the available phenomenological representations 
of the scattering amplitude. 



12 Subthreshold expansion 

Hohler and collaborators |2^ analyze the low energy structure in terms of 
an expansion of the amplitude around the point u = t = 0. More precisely, 
the expansion concerns the difference between the full amplitude and the 
pseudovector Born term, specified in eq. Crossing symmetry implies 

that the amplitudes X G {D~^ , D~ /u, /u, B^} are even in i/, so that the 
expansion takes the form 

X{iy, t) = Xp^{iy, t) + Xqq + Xio + Xqi t + X20 + Xxv v^t + X02 + • • • 
At leading order, only the tree graphs from contribute: 

4o = 0{M^) , doo = ^ + 0{M^) , 600 = ^ + 0(1) . 

The first two relations may be viewed as a variant of Weinberg's predictions 
for the ^-wave scattering lengths - the two low energy theorems discussed 
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in section |Tl] represent yet another version of these relations. The third 
does not contain significant information, because the coefficient fogo picks up 
contributions of the same order also from C}^ - chiral symmetry does predict 
the values of d^^ and c/qq for = = 0, but does not constrain 6qq. 

The representation of the scattering amplitude to O^q"^) is fully deter- 
mined by the coupling constants of the effective theory. Accordingly, we can 
calculate the various coefficients in terms of these, to the corresponding order 
of the chiral expansion. It suffices to expand the loop integrals first around 
ly = t = and then in powers of M. For c/qq, for instance, the calculation 
yields 



^00 



2F3 



1 + Ml {8id, + d2 + 2d,) + ^[X^ + 



3F3 \ 32n 



{ 64:71 F^ttin 2 7rF^ J 

We have expressed the result in terms of the physical values of F^^, M^^ and 
Qa. The formula represents an exact result: Up to and including terms 

3 

of 0(m2), the expansion of d^Q in powers of the quark masses is fixed by 
the coupling constants of the effective Lagrangian specified in section ^. In 
particular, the result shows that this expansion contains a chiral logarithm 
at order M^, as well as a term of O(M^). Similar formulae can be given also 
for the other subthreshold coefficients. These are listed in appendix ^ Note 
that, for some of the coefficients, the expansion starts with M~^, indicating 
that those coefficients diverge in the chiral limit, on account of the infrared 
singularities required by unitarity. 

In the loop integrals, the coupling constants ci, ... , C4 are needed only to 
leading order. As discussed in detail in ref. we may invert the relations 
for (ioo, d'^Q, c/qi, ^qq and represent these couplings in terms of the subthreshold 
coefficients: 

ci = (4o + 2 Ml + 0(M.) , c, = ^dto + 0(M.) , (12.1) 

C3 = -Fl < + 0(MO , C4 = 600 - + 0{M^) . 

Apart from the Born term, the tree graph contributions to the invariant 
amplitudes represent polynomials in the momenta. The coefficients of these 
polynomials involve precisely the same combinations of effective coupling 
constants as the corresponding subthreshold coefficients - we may trade one 
set of parameters for the other: Up to and including 0(g*^), the scattering 
amplitude may equally well be expressed in terms of subthreshold coefficients. 
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13 Chiral symmetry 



When parametrizing the scattering amphtude in terms of the subthreshold 
coefficients, part of the information is lost: Chiral symmetry imposes con- 
straints on these coefficients. In order to identify the constraints, it is useful 
to split the effective Lagrangian into a chirally symmetric part, Cg-, and a 
remainder, Csb-, that collects the symmetry breaking vertices. The leading 
term C^^^ does not contain a symmetry breaking piece. In C^m\ the term pro- 
portional to Ci(x+) belongs to Cgh, while those involving 02,03,04 survive in 
the chiral limit and hence belong to At order g^, the coupling constants 
(is, diQ, dis represent symmetry breaking effects. At 0{q'^), we did not specify 
the Lagrangian, but instead wrote down the corresponding contributions to 
the scattering amplitude, in eq. ( [4.11 ). Concerning D~^, only those terms in 
that contain four derivatives of the meson field belong to the symmetric 
part. In view of 2g ■ g' = 2M^ — t and Amj^i' = [P + P') ■ [q + q'), the 
expressions (t — 2M^)^, {t — 2M^) z/^ and z/^ are of this type. Hence the 
coupling constants cq, 67, cs, belong to the chirally symmetric part of the 
Lagrangian, while the combinations 63 — 4e8, 64 + 2e7 and 65 + Aeg represent 
symmetry breaking terms. For the same reason, cio, Cn belong to Cs, while 
the combination eg + 2eii controls the symmetry breaking in the amplitude 

. The coupling constants Ci and 62, that occur in the chiral expansion of 
the scalar form factor at ffist non-leading order, also belong to the symmetry 
breaking part of the Lagrangian. 

If the quark masses are turned off, only the coupling constants of Cg con- 
tribute. One readily checks that, in that limit, the subthreshold coefficients 
listed in appendix ^ represent linearly independent combinations of the sym- 
metric coupling constants, except for two constraints: d^Q vanishes and d^Q is 
determined by F - precisely the two low energy theorems encountered already 
at leading order. We conclude that, in addition to the Goldberger-Treiman 
relation, chiral symmetry only imposes two conditions on the subthreshold 
coefficients. 

As noted above, we can express the amplitude in terms of subthreshold 
coefficients instead of coupling constants, using the relations (|12.1|) to elim- 



inate Ci, . . . , C4 in the contributions from the loops. For the value of F^D^ 
at the Cheng-Dashen point, we then obtain 

S = (k, rf+ + 2 K^M^d^, + AMtdt,, + K^Mldt^) + A 



4o + 2 K^Mld^^ + 4M^42 + K^Mldt^) 
1 q A/f2 9. A/f2 



K2 = 1 + 



96 7r2F2 64 ttF^ 
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128 7rF2 ' 



960 7r2F2 

3840 vrFJ mn'^m^F^ ' USixm^F^ 



+ 



+ 0{Ml] 



The only infrared singularity occurring in this relation is the term propor- 
tional to M^. Since the coefficient of this term is very small, the contributions 
from the loop graphs are tiny: Numerically, we find Ki = 1.003, K2 = 1.012, 
K3 = —0.001 and A"*" = 1.1 MeV. Inserting the values of the subthreshold co- 
efficients from ref. the relation yields E = 61 MeV, in perfect agreement 
with the value given in that reference. 

We conclude that, in the present case, the expansion is rapidly convergent. 



According to section 11, the term E nearly coincides with the value of the 
scalar form factor at t = 2M^. The difference between a{2M^) and a = a{0) 
is well understood - the evaluation within chiral perturbation theory |^ 
confirms the result of the dispersive calculation described in ref. ||23|| : 



(t(2M^) - (t(0) = 15.2 ± 0.4 MeV 



'13.21 



The relation ( |13.1| ) may thus be rewritten in the form 

a = Fl (4o + 2 Midi, + AMtdt^) - . (13.3) 
Evaluating the small correction terms with the subthreshold coefficients of 



ref. f2Sl, we obtain 



di = 13 ± 2 MeV , 



;i3.4) 



where the error is our estimate for the uncertainties arising from the correc- 
tion Aci3, as well as from the higher order contibutions, which our calculation 
neglects. 

For the isospin odd amplitude D~ , the term analogous to E is the quantity 
C, defined in eq. ( 11. 5| ). The representation in terms of the subthreshold 
coefficients is of the form 



C = 2F2(rf- +2MX-i) + A-, 
where the correction A~ is given by 

{hgl-2)Ml {gl-l)Ml gl{3g 



(13.5) 



22) M3 



+ 



72 7r2F2 



lQV2 7rm^F^ 



96 7rF2 192 7rm^F2 

ln(l + V2) + 0{Mt) . 
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The infrared singular terms proportional to are small, of order 0.03. 
Moreover, they nearly cancel against those of 0{M^), so that the net cor- 
rection from the one loop graphs is tiny, A~ = —0.003 + O(M^). The con- 
tributions of order are more important. In particular, the A resonance 
generates a term of this order, which is enhanced by the third power of the 
corresponding energy denominator (note that the resonance also contributes 
to the coefficients d^Q and (Iqi): 

, _ _ SF^gjM:: 
^ 9mjv(mA — mjv)^ 



Numerically, this term amounts to 0.02.0 The dispersive analysis of ref. ||28 
confirms that the resonance contribution dominates. In our opinion, the 
estimate 

A- =0.02 ±0.01 (13.6) 

covers the uncertainties. We conclude that the relation ( |13.5| ) very accurately 
determines the value of the isospin odd amplitude at the Cheng-Dashen point, 
in terms of the subthreshold coefficients c^qq, d^i- 

Unfortunately, our knowledge of the coefficient rfog is not satisfactory. 
The result obtained for this quantity is very sensitive to the value used for 
g-^N- The difference d^Q — glj^/2m^ is relatively well known, as it is given by 
an integral over a total cross section 



^oo-7r-T = - / -^ka (ij). (13.7) 



If the integral is taken known, a shift in = gl^M^/lGnm'^ from the value 
P = 0.079 [gD to /2 = 0.076 O lowers the result for cioo by O.OSikf-^, so 
that C decreases by 0.07. Indeed, the recent evaluations of the subthreshold 
coefficients by Stahov [^| and Oades confirm that the value of c/qq is 



lower for the partial wave analyses of the VPI/GW group than for the Karl- 
sruhe analyses. For a precise determination of the value of the isospin odd 
amplitude at the Cheng-Dashen point, the discrepancy in the value of the 
pion nucleon coupling constant needs to be resolved. 

We can now check whether the phenomenological information is consis- 
tent with the low energy theorem ( |11.6|) , which predicts the value of C to 
order M^, except for the contribution Cgb = ^^M^dKfi) from the symme- 
try breaking coupling constant d^. For the reasons given earlier, we expect 
contributions of this type to be very small, of the order of a few per cent. 

^ The corresponding contribution to A+ is even smaller, of order 0.1 MeV. 
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^00 




C 




KH80 [H 


1.53 ± 0.02 


-0.134 ± 0.005 


1.13 ± 0.02 


0.04 ± 0.03 


KA84 g 


1.510 ± 0.001 


-0.136 ± 0.003 


1.11 ± 0.01 


0.03 ± 0.01 


SP98 g 


1.468 ± 0.004 


-0.138 ± 0.003 


1.06 ± 0.01 


0.00 ± 0.01 


KH80 [0 


1.53 ± 0.02 


-0.14 ± 0.03 


1.12 ± 0.06 


0.04 ± 0.06 


SP98 g 


1.46 ± 0.01 


-0.12 ± 0.01 


1.09 ± 0.02 


0.06 ± 0.02 



Table 1: Value of the isospin-odd amplitude at the Cheng-Dashen point. The 
results for C are obtained with eqs. (|13.5| ) and ( |13.6| ). The value of d^fi) 



follows from eq. (|11.6|) and refers to the running scale fi = 1 GeV. 



This is confirmed by the numerical results listed in the last column of table Q 
which are obtained by comparing eqs. ( |11.6| ) and ( |13.5| ) - the relation ( |12.1| ) 



is used to express Ci in terms of subthreshold coefficients, the running scale 
fi is set equal to 1 GeV and contributions of order or higher are dropped. 

As mentioned above, the Weinberg predictions for the two S-wave scat- 
tering lengths represent yet another version of the two energy theorems under 
discussion. These predictions only hold to leading order of the chiral expan- 
sion, but the one loop representation of the amplitude allows us to analyze 
the corresponding corrections - it suffices to evaluate our representation of 
the amplitudes D^{v,t) at u = M^^, t = 0. The expansion of the result in 
powers of then yields an explicit expression for the S-wave scattering 
lengths, valid up to and including O(M^). In fact, this calculation has been 
done already, in the framework of HBCHPT and we have verified that 
our method yields the same result (note that the results given in Appendix 
A of 1^ contain typographical errors). 

In the case of the scattering lengths, however, the corrections are by no 
means small - the first few terms of the chiral expansion only yield a rather 
poor approximation. This is in marked contrast to the case of vrvr scattering, 
where the chiral expansion of the scattering lengths converges rather rapidly. 
The main reason is that the distance between the Cheng-Dashen point and 
threshold is of order M^^ and is therefore much larger than for vrvr scattering, 
where the analog of the Cheng-Dashen point occurs at s = u = M^, t = 
2M^, so that the distance to threshold is of order M^. Accordingly, the 
expansion parameter is larger - we cannot expect the chiral representation 
at threshold to have the same accuracy as in the case of vrvr scattering. The 
problem is accentuated by the fact that the threshold sits on top of the 
branch point required by unitarity - when considering the scattering lengths, 
we are in effect analyzing the amplitude at a singular point. In this respect. 
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the situation is comparable with the case of tttt scattering, where the 5- 
wave scattering lengths also contain chiral logarithms with a large coefficient. 
As shown in ref. |T^, the expansion of the tttt scattering amplitude in the 



subthreshold region converges much more rapidly than at threshold. 

In the variants of the low energy theorems discussed above, these prob- 
lems did not occur, because we have been comparing the properties of the 
amplitude at the Cheng-D ashen point, s = u = m^, t = 2M^, with those 
at s = u = + Ml, t = 0. In this region, the amplitude does not contain 
branch points and the relevant distance is small, of order M^. This is why 
the corrections encountered turned out to be rather small. 



14 Cuts required by unitarity 




-4-2 2 4 



Figure 7: Mandelstam plane. Hatched: Physical regions of vrA^ scattering. 
In the shaded region s < {rrii^ + M^)^, u < {rrii^ + M^r)^ and t < AMI the 
amplitude is real and analytic, once the poles at s = and u = (thin 
lines) are subtracted. Dotted: Nucleon poles. Dashed: A-resonance. Dot- 
dashed: p-resonance 
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In the following sections, we examine the energy dependence of the loop 
graphs in some detail, in order to clarify the problems encountered when 
extrapolating the chiral representation to threshold or even beyond. For this 
purpose, we first study the analytic structure of the loop graphs. 

In addition to the nucleon pole terms, the scattering amplitude contains 
branch points at s = (m^r + M^)^, u = {m^ + M^r)^ and t = 4M^. These 
correspond to the sides of the shaded triangle shown in fig. 0. In the chiral 
representation of the amplitude, the corresponding cuts are described by the 
loop integrals. Crossing symmetry relates the discontinuities in the w-channel 
to those in the s-channel, so that we only need to discuss the s- and t-channel 
imaginary parts. 

We recall that we are absorbing the self energy insertions due to the 
tree graphs of C^^^ and by replacing the bare nucleon mass with m4 = 
m — 4ciM^ + ei M^. As far as the loop integrals are concerned, the difference 
between 777.4, M and the physical masses rriff, M^, is beyond the accuracy of 
our calculation. In the following, we represent these integrals in terms of the 
physical masses. 

14.1 s-channel cuts 

The branch cut at s > (777,jv + M,,y is an immediate consequence of the 
unitarity condition, 

{f\T\t) - (/|Tt|^) = ^5^(/|Tt|77)(77|T|/) . (14.1) 

{n} 

At one loop level only nN intermediate states enter this relation - inelas- 
tic channels only show up at higher orders. The condition thus exclusively 
involves the s-channel partial wave amplitudes 

fLis) = q-'expt6^^is)sm6l^is) . (14.2) 

The lower index specifies the orbital angular momentum I as well as the total 
angular momentum J = l±^. The upper index refers to the isospin = |, f 
and q stands for the cm. momentum, 

g2 = 1 {s - (777^ + M^y}{s - {m^ - M^f} . (14.3) 
4s 

In this normalization, elastic unitarity takes the form 

Im/4(s)=g|/4(s)r (14.4) 
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The low energy expansion of the partial waves starts at 0{q). The relation 
(|14.4| ) implies that an imaginary part starts showing up only at 0{q^) - the 
tree graphs generated by and thus suffice to evaluate the imaginary 
parts of the amplitude to 0{q^). In agreement with the chiral counting rules, 

(2) 

the imaginary parts generated by the square of the contributions due to C\, 
are of 0{q^) and we can drop these. 

We have checked that our representation of the scattering amplitude in- 
deed obeys elastic unitarity. In this respect, there is no difference between 
infrared and dimensional regularization - the s-channel discontinuities are 
identical 0. The real parts are different, however. To any order in the 
low energy expansion, the difference is a polynomial, but the degree of the 
polynomial increases with the order. 

14.2 t-channel cuts 

Some of the graphs contain a branch cut for t > 4M^. The corresponding 
discontinuity is determined by the extended unitarity relation 

Im fiit) = {1 - 4M2/t}^ t'j{ty fiit) , 4M2 < t < iml . (14.5) 

The quantities f±{t) represent the t-channel partial waves of the vrA^ scatter- 
ing amplitude. The quantum number J stands for the total angular momen- 
tum, while the lower index again refers to the spin configuration. There is 
no need to in addition specify the isospin quantum number, because it is de- 
termined by the total angular momentum: J even —>■ I = 0, J odd ^ / = 1. 
The functions 

t'jit) = {l-AM'Jt}-hxp^6'Jit) sm5'jit) 

denote the partial wave projections of the tttt scattering amplitude. 

The low energy expansion of the right hand side of ( |14.5| ) again starts at 
O(g^). In our context, the evaluation of this condition is particularly simple, 
because the vrvr scattering amplitude is needed only to leading order, where 
it is a polynomial. Hence only the partial waves with J = 0, 1 contribute. 
These are given explicitly in appendix 0. 

We have checked that, in the region t < 4m^, the imaginary parts of our 
loop integrals do obey the relation ( |14.5|) . Note that this relation only ac- 
counts for the TTTT intermediate states. States with three or more pions do not 
occur at one loop, but outside the low energy region, the nucleon-antinucleon 
states yield an additional contribution, which does show up in the box graph 
and in the triangle graph, fig. Qm, through a branch cut for t > 4m^. Be- 
cause this singularity is absent at any fixed order of the chiral expansion of 
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the loop integrals, it is not preserved by infrared regularization. Indeed one 
finds that the regular part of the corresponding loop integrals has a branch 
cut for t > 4m^ so that the imaginary part of the amplitudes in infrared reg- 
ularization agrees with the one obtained in dimensional regularization only 
for t < imt. 



14.3 Box graph 

The analytic structure of the Box graph (i) is more complicated: The corre- 
sponding loop integrals contain a simultaneous cut in the variables s and t 
and thus contribute to the Mandelstam double spectral function (at one loop 
level, only the box graph generates such a contribution). 

Let us consider the relevant scalar loop integral. In view of the four 
energy denominators, this integral does not require regularization. It may 
be represented in terms of the double dispersion relation 

The boundary of the integration region is given by 

[s - s+) [s - s_) 
and the explicit expression for the double spectral density readsQ 



8a/ t (t - C,„(s)) (s - s+) (s - s_) 

In view of t„,i„(s) > 4m^, we may expand the integral in a Taylor series 
in powers of t. Because the singularity in the variable t lies far outside 
the low energy region, the higher order terms in t are strongly suppressed 
- these contributions are similar to those we dropped when simplifying the 
integrands of the loop integrals in section ^ The statement also holds in 
infrared regularization and for the tensor integrals associated with the box 
graph. 

Using the explicit expressions for the box graph in appendix 0, we have 
determined the order to which the t-dependence of the integrals /i3(s,t), 

^The expression for p2{s,t) may be read off from the one given for the imaginary part 
Ims Hi3{s, t) = luis Iisis, t) in appendix The square root occurring in that expression 
becomes imaginary for t > tmi„- The double spectral density is the discontinuity across 
the cut generated by this square root. 
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l[]^{s,t), l[f{s,t) is needed - because of cancellations between the various 
terms occurring in the representation of the invariant amplitudes in terms of 
these, the order needed cannot directly be read off. Within the accuracy of 
our calculation, the box-amplitudes of are given by their value at t = 
and the amplitudes can be replaced by their first order expansion in t, 
so that we are allowed to replace the box integrals by their Taylor expansion 
in t to first order. 

The expanded box integrals can be represented in terms of a spectral 
function with a single variable. In the case of the full scalar integral, for 
instance, the explicit representation reads: 

TT s' - s-ie ^ 3C„(s) ^ '*) 

v/(g-g-)(g-g+) 

167r(s -ml- 2M2) {mis - {ml - Mlf) ' 

The formula nicely displays the suppression of the t-dependent part. Since 
tmin{s) is a quantity of order g° it is also evident that the t-dependent part 
is an 0{q^) correction to ifi3(s,0). 

This shows that, at low energies, the contributions generated by the Man- 
delstam double spectral function can be described in terms of single variable 
dispersion integrals. 

15 Dispersive representation 

With the above simplification of the box graph contribution and after the 
reduction of the numerators for the diagrams discussed in section ^ all of 
the loop integrals are of a similar structure: The discontinuities across the 
s- and M-channel cuts are linear in t, while those across the t-channel cut are 
linear in v. Up to higher order corrections, the s-channel imaginary parts 
thus have the form 

Im,D^(s,t) = lmDf{s)+tlmD^{s) + 0{q^) , 
Im,5±(s,t) = lmBf{s) + 0{q^) , 

while the t-channel discontinuities have the structure (explicit expressions 
for the functions ImDj(t), ImD^(t), Imi?^(t) are given in appendix 0) 

lmtD+{v,t) = lmD+{t) + 0{q^) , lmtD~{u,t) = ulmD^ {t) + 0{q^) , 
ImtS+(z/,t) = 0{q^) , lmtB-{iy,t) = lmB^{t) + 0{q^) . 



Hi3{s,t) = 
Pi(s) = 
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We now make use of the fact that an analytic function that asymptot- 
ically only grows with a power of the argument is fully determined by its 
singularities, up to a polynomial. This allows us to establish an alternative 
representation of the scattering amplitude that also holds to first nonleading 
order and has two advantages as compared to the one in terms of the explicit 
expression for the loop integrals: It exhibits the structure of the amplitude in 
a more transparent manner and only involves the expressions for the imag- 
inary parts. Since unitarity determines the latter in terms of the tree level 
amplitude, the dispersive representation does not involve any loop integrals. 

Explicitly, the dispersive representation reads 

D+{u,t) = D;^{u,t) + Dtis) + D+iu) + tDt{s) + tDt{u) 

+ DUt) + D;iu,t) + Oiq'), 

D-{u,t) = D-^iu,t) + D^{s)-D^{u)+tD,{s)-tD,{u) (15.1) 

+ uD^{t) + D^iu,t) + 0{q') , 

B+{u,t) = B+{u,t) + B+{s)-Bt{u) + B^{u,t) + 0{q^) , 

B-{u,t) = B;^{u,t) + B^is) + B^iu) + B,{t) + B;iu,t) + 0{q^) . 

The functions D^^{u,t), B^^iy^t) represent the pseudovector Born terms of 
eq. The contributions generated by the cuts in the s-, t- and ^-channels 

are described by the functions D^^s), Df{u), . . . , -82^ (t). The dispersion in- 
tegrals that determine these in terms of the corresponding imaginary parts 
are given below. Once the singularities are removed, the low energy repre- 
sentation of the invariant amplitudes reduces to a polynomial. We denote 
the corresponding contributions by (z/, t), _B^(i/, t). 

A representation of the above type was used already by Gasser, Sainio 
and Svarc [0, |3^. These authors noted that if the masses occurring in the 
loop integrals are identified with the physical masses rather than the bare 
quantities occurring in dimensional regularization, the infrared singularities 
generated by the cuts are correctly accounted for. As they did not make 
use of any approximations for the loop integrals, the t-dependence of their 
s- and M-channel cuts was however more involved. The representation (|15.1|) 
shows that, to the accuracy of a one loop calculation, the amplitude may 
be described in terms of nine functions of a single variable. An analogous 
representation also holds for the tttt scattering amphtude In that case 
three functions of a single variable suffice and the representation is valid even 
to two loops, that is modulo contributions of O(g^). 
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16 Subtracted dispersion integrals 



We now specify the functions Df(s), . . . ,Bf(t) as suitably subtracted dis- 
persion integrals over their imaginary parts. For this purpose, it is convenient 
to replace the variable s by the lab energy uj, which is linear in s and inde- 
pendent of t (for t = 0, u coincides with u), 

. = (16.1) 

2miv 

The number of subtractions is determined by the asymptotic behavior of the 
imaginary parts and is related to the superficial degree of divergence oj{G) 
of the relevant diagram. This quantity is obtained by adding up the powers 
of the loop momentum k in the region where it is large compared to the 
masses and to the external momenta: oj{G) = 4 + 6 — Ij^ — 2 1^^, where 6 is 
the number of derivatives at the vertices and Jjv, J^r stand for the number 
of nucleon and pion propagators occurring in the loop, respectively. The 
vertices from contain one derivative, those of C^'^ contain two. In C^^\ 
there are up to four derivatives (vertex proportional to C2). It is important 
here that we are considering only diagrams with at most one vertex from 
C\f - if the unitarity relation is evaluated with the full square of the tree 
level amplitude, the imaginary part grows more rapidly. Note also that the 
simplifications discussed in section ^ reduce the degree of divergence. 

The imaginary parts occurring in the tensorial decomposition of the loop 
integrals are listed explicitly in appendix |^. These expressions show that 
the quantities ImDf(s) grow with s'^lns, while for the parts linear in t, we 
have ImDf(s) oc sins. The leading contribution stems from the graphs 
(g) and (h) and is proportional to C2 (if only the loops generated by 
are considered, we instead obtain lmDf{s) oc s^ and ImDf'(s) oc s). The 
asymptotic behavior of the ^-amplitudes is given by Imi?^(s) oc C2slns, 
Imi?f(s) oc s. We subtract at = and write the s-channel dispersion 
relations in the form 

, p duj'lmDfis 

(s) — — 



Jm^ uj'"^{uj' -u- ie) 



Bt{s) -I 



vr ./m, ^'^ {uj' -iu- ie) 



The representation ( |15.1 ) is manifestly crossing symmetric: The contribu- 
tions generated by the cuts in the w-channel are are obtained from those due 
to the s-channel discontinuities by the substitution s — > m. 
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The contributions due to the t-channel cuts may be characterized in the 
same manner. The dispersion relations take the form 



p dt'lmD+{t') 



TT 74JV/2 t'^ {f — t — it) 



e dt' ImD7(t') 



t p dt'lmB^it') 



TT JiM^ t' {f — t — it) 



In the terminology used in ref. the contributions associated with the 
cuts in the s-, t- and w-channels represent "unitarity corrections". Together 
with the explicit expressions for the imaginary parts of the various graphs 
given in the appendix, these relations unambiguously specify the contribu- 
tions generated by the cuts. 

By construction, the cuts do not contribute to the first few terms in the 
Taylor series expansion of the amplitude in powers of v and t, that is to 
the leading coefficients of the subthreshold expansion, which we discussed in 
section 12. The polynomials i?^(z/, t) and Bp[u,t) account for the following 
terms of this expansion: 

D+{v, t) = d^Q + rf+ z/2 + rf+ t + d^Q + rf+ vH + rf+ , 

B;{v,t) = btoU, (16.4) 

Dp{u, t) = d^f^u + d^Q + c?oi i/t , 

Bp {u, t) = 600 + bio '^^ + boit ■ 

At the order of the chiral expansion under consideration, the higher coeffi- 
cients exclusively receive contributions from the parts of the amplitude that 
describe the cuts. 

To complete the construction of the dispersive representation, we incor- 
porate the constraints imposed by chiral symmetry. These are: 

1. The strength of the vrA^ interaction is determined by g^, up to correc- 
tions from the symmetry breaking part of the effective Lagrangian. 

2. In view of the relations ( |12.1| ), chiral symmetry and unitarity fix the 
imaginary parts in terms of the subthreshold coefficients. 

3. Finally, as shown in section O, chiral symmetry subjects the subthresh- 
old coefficients to the two constraints ( |13.1|) and (|13.5|) . 
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With these supplements, the dispersive representation of the amphtude be- 
comes equivalent to the one in terms of infrared regularized loop integrals. In 
effect, we made use of these integrals only when analyzing the corrections to 
the Goldberger-Treiman relation and to the two low energy theorems referred 
to in point 3. 

Note that the dispersive representation does not rely on an expansion of 
the infrared singularities in powers of the momenta - as mentioned above, 
this is a subtle matter. In particular, the result depends on the choice of the 
two independent kinematic variables that are kept fixed when performing the 
expansion. The above definition of the cut contributions in terms of their 
absorptive parts precisely serves the purpose of avoiding such problems at 
this stage. We will give the explicit algebraic representations for the functions 
D^{s), . . . , Bf{t) that result, for instance, if the expansion is performed by 
keeping the ratio oo/M^^ fixed and discuss the changes occurring if we replace 
the lab energy u by the cm. energy. Also, we will discuss the range of validity 
of such representations in some detail, but first we wish to describe the result 
of our calculation without invoking expansions of this sort. 

17 Dispersion relations versus loop integrals 

The dispersive representation of the loop integrals may be viewed as yet an- 
other regularization of the loop graphs that is consistent with power counting. 
We add a few remarks concerning the comparison of this representation with 
the one based on infrared regularized loop integrals. For this purpose, we 
first again consider the scalar self energy integral, for which both of these 
representations can be given explicitly (compare section |5.1|) . In dimensional 
regularization, the relevant integral is given by H{s), while in infrared regu- 
larization, only the infrared singular part thereof, I{s), is retained. On the 
other hand, the dispersive representation of the self energy integral, which we 
denote by D{s), is determined by the imaginary part of I{s). The relevant 
dispersion relation involves the minimal number of subtractions required by 
the asymptotic behaviour: 



On the interval of interest, s > (m^v+MTr)^, we have ImJ(s) = Imif(s). While 
the function H{s) is analytic except for the discontinuity on this interval, the 
function I{s) in addition contains fictitious singularities far outside the low 
energy region: a pole at s = 0, as well as a cut along the negative s-axis. 
Since by construction, the dispersive representation only contains the right 
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hand cut, D{s) can differ from the function H{s) only by a polynomial. 
In fact, the asymptotic behaviour shows that the polynomial is a constant. 
Hence the explicit expression for D{s) reads 



D{s) = H{s) - H{ml + M^^ 

In view of H{s) = I{s) + R{s), the difference between the two representations 
we wish to compare may be represented as 

Dis) - lis) = Ris) - R{ml + Ml) - I{ml + Ml) . 

Since the regular part of the self energy integral, -R(s), does not contain any 
infrared singularities, we may replace it by the first few terms in the chiral 
expansion, which in the present case amounts to a simultaneous expansion in 
powers of uj and M^. The resulting representation for R{s) is a polynomial 
in UJ, with coefficients that are polynomials in Ml. We conclude that, to 
any finite order in the chiral expansion, D{s) — I{s) is a polynomial in uj or, 
equivalently, in s. Moreover, since all of its coefficients except the constant 
term only pick up contributions from -R(s), they are polynomials in M^. The 
constant term, on the other hand, does contain infrared singularities, from 
I{ml + Ml). 

The same line of reasoning also applies to the other graphs. When form- 
ing the sum of all contributions, the difference also amounts to a polynomial 
and only the coefficients listed in eq. ( |16.4D can pick up infrared singularities. 



from the analogues of the term J(m^ + M^) encountered above. As we are ex- 
plicitly supplementing the dispersion integrals with the relevant polynomial, 
the result agrees with the representation in terms of infrared regularized loop 
integrals, up to terms that are beyond the accuracy of a one loop calculation. 



18 Chiral expansion of the amplitude 

The dispersive representation given above includes terms of arbitrarily high 
order in the chiral expansion. It is suitable for numerical analysis, but since 
the expressions for the imaginary parts are lengthy, an approximate represen- 
tation that replaces these integrals by elementary functions is useful. Such 
a representation may be obtained by performing the chiral expansion in the 
integrand of the loop integrals. The result is very similar to the one obtained 
in HBCHPT, where the unitarity cuts are described in terms of simple func- 
tions such as arccos(x). In 0, we have proven that this expansion converges 
for all of the integrals encountered here, at least in part of the low energy 
region. In the present section, we first briefly comment on the chiral expan- 
sion of the integrals associated with the t-channel cut, where the discussion 
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given in ref. can be taken over without significant complications. Then 
we analyze the analogous issue for the s-channel and discuss the infrared 
singularities associated with the corresponding cut in some detail. 



18.1 t-channel cut 

Only the diagrams belonging to the topologies (i), (k), (1) and (m) develop 




(i) (k) (1) (m) 

an imaginary part for t > 4M^. For graphs (i) and (m), the t-channel cut 
only starts at t = 4m^. As discussed in section |T^, these do not show up at 
finite orders of the chiral expansion, so that the diagrams (i) and (m) do not 
contribute. The integrals arising from the graphs (k) and (1) are polynomials 
in u. It is convenient to replace the pion mass and the momentum transfer 
by the two dimensionless variables 

_ _ t 

Graph (k) may be expressed in terms of the scalar loop integral with two 
meson propagators, J(t). We remove the divergence by subtracting at t = 0. 
The remainder, J(t) = J{t) — J(0), only depends on the ratio r and the 
explicit expression reads 



Jir) = — - < 1 — -i / arcsm — 

^ ^ 87r2 I V r 2 

The triangle graph (m) involves integrals with two meson propagators and 
one nucleon propagator. The chiral expansion of these functions does not 
cover the entire low energy region: It breaks down in the vicinity of the 
point t = 4M^. The matter is discussed in detail in ref. ||^, where it is shown 
that the elementary function 



9(r) = ^ln^-^ln|l " 



32 7iy/T 2-v^ 32 TT t V4^ 




327r2 1^ ^/4^ v/r(4-r) 
does provide a representation that is valid throughout the low energy region. 
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18.2 s-channel cut 



The dispersive integrals associated with the s-channel cut are polynomials in 
t with coefficients that depend only on s. The w-channel cut is described by 
the same integrals. The effective theory concerns the region where s — = 
0{q). We may for instance consider the expansion of the loop integrals in 
powers of a = MT,/mj^ at a fixed value of the ratio r = (s — m^)/M^. That 
expansion, however, necessarily diverges in the vicinity of the branch point 
s = {nif^ + M^)^, because this singularity corresponds to a value of r that 
depends on the expansion parameter a. A better choice of the variable to be 
held fixed is the ratio of the pion lab. energy uj to the pion mass, which we 
denote by 

^ s — ml — Ml UJ 



2 

The s-channel cut starts at = 1 and thus stays put when a varies. Indeed, 
we have shown |P that the convergence domain of the expansion of the self 
energy integral I{s) at fixed is < s < 2m^ + 2M^ and thus covers exactly 
the Mandelstam triangle. 

Numerically, the convergence is rather slow, however, already at the 
threshold. The explicit expression for the self energy integral I{s) involves 
the factor m? / s = {1 + 2aQ + and some of the loop integrals even 

contain the third power thereof. At threshold, the expansion of this fac- 
tor in powers of a is roughly a geometric series with expansion parameter 
2Qa ~ |. When the factor appears squared or cubed, the lowest orders do 
therefore not give a decent approximation. 

A more suitable choice is the c. m. energy of the pion 

= = ^ , (18.1) 

which maps the interval < s < oo onto — oo < < oo and thus pushes 
the singularity at s = away from the threshold at Qg = 1. Expressed in 
terms of this variable, the infrared singular part of the self energy integral is 
given by 

I{s) = {- ~ arccos(-l]g) + 1^, (l - 32 tt^ A^) | . (18.2) 



At fixed Qq, only the factor ^/s = y 1 + (Q^ — l) + aQq gets expanded. 

It has a branch point at Q'^ = — ja. The expansion is convergent 

for < which corresponds to the range 

(3 - 2^){ml - Ml) < s < (3 + 2^){ml - Ml) . 
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On the left, the convergence range is shghtly smaller than the one obtained 
for the variable Q, but on the right, the convergence is drastically improved. 
The point at which the expansion at fixed fig breaks down is six times farther 
away from threshold than the edge of the Mandelstam triangle, where the 
expansion at fixed Q starts to fail. 




Figure 8: Chiral expansion of the contribution from diagram (f) to the real 
part of the amplitude D'^ (infrared regularization, scale fi = itim)- Full: 
relativistic amplitude, dot-dashed: expansion at fixed fig, dashed: expansion 
at fixed VL. In all of the figures, dimensionful quantities are given in units of 

We illustrate this with the contribution of graph (f) to the amplitude 
at t = 0. The corresponding Feynman diagram is shown in fig. |^ and the for- 
mulae in appendix explicitly specify this contribution in terms of the func- 
tion I{s). In fig. ^, the full line shows the exact result for the renormalized 
amplitude at t = 0, obtained with infrared regularization (for definiteness, 
we have identified the running scale with the nucleon mass). The dashed line 
corresponds to the expansion to order a'^ at fixed f2, while the dash-dotted 
curve is obtained if we instead keep the variable VLq fixed. The figure shows 
that the chiral expansion at fixed Vtq describes the result quite accurately. 

In fig. 1^, we compare the imaginary parts of the loop integrals with the 
chiral expansion thereof, for the sum of all graphs. Near threshold, these are 
suppressed by phase space. To make the behavior visible there, we plot the 
integrand relevant if the variable in the dispersion integrals is identified with 
the lab momentum k = ^Juj"^ — rather than the lab energy u. The inte- 
grand of the dispersion relation for the value of the amplitude at threshold, 
D^{M^, 0) = 4 7r (1 + a) a^+, is then given by n^{k) lmDf{k), with 



TTUJ^ dk 1 - uj + \ Trkuj'^ [ (jJ 
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Figure 9: Imaginary part of the D-amplitudes at t = 0. The normahzation 
factor n^{k) is chosen such that the area under the curve represents the 
contribution to the value of the amphtude at threshold. Full: relativistic 
amplitude, dot-dashed: expansion at fixed fig, dashed: expansion at fixed Vt. 



In view of the difference in the weight, the dispersion relation for D is more 
sensitive to the high energy region than the one for . 

18.3 Chiral expansion of the dispersion integrals 

The situation changes drastically, however, if we consider the correspond- 
ing contribution to the subtracted dispersion integral Di{s), which may be 
evaluated as follows. We first note that the discontinuity across the right 
hand cut of the function I{s) coincides with the one of the corresponding 
dimensionally regularized object, H{s). Moreover, H{s) does not have a left 
hand cut. Hence the expression obtained by simply replacing I{s) with H{s) 
does have the same singularities as the dispersive integral we are looking 
for - the difference is a polynomial in u. The subtracted dispersion inte- 
gral is obtained from H{s) by removing the Taylor series expansion to order 
cj^. In fig. 0, the full line represents the exact contribution to IieDf{s), 
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while the dot-dashed one depicts the chiral expansion thereof at fixed fig, 
to 0{a^). Evidently, the truncation of the chiral expansion at O(a^) very 
strongly distorts the result. 



ReD+(s) 




Figure 10: Chiral expansion of the contribution from diagram (f) to the 
subtracted dispersion integral D^{s). 



The reason for this disaster is that the chiral expansion of the Taylor co- 
efficients that specify the subtraction polynomial converges only extremely 
slowly. The contributions from the s-channel loop graphs to the coefficient 
d2o, for instance, are obtained by evaluating the fourth derivative of the cor- 
responding contributions to at u; = 0. These quantities are represented 
very poorly by the leading terms in their chiral expansion. The problem 
disappears if we keep the Taylor coefficients fixed when performing the chi- 
ral expansion (the representation that then results is indicated by the dotted 
line). When representing the amplitude in the form ( |15.1D , we have absorbed 
the subtraction polynomial in the coefficients of the subthreshold expansion. 
This implies that the chiral expansion of the subthreshold coefficients in pow- 
ers of the quark masses - in particular the one of d^g ^ converges only very 
slowly. 

The physics underneath the phenomenon is readily identified. The coef- 
ficient '^^y tie expressed as an integral over the total cross section. The 
integral is dominated by the threshold region. If the quark masses are sent 
to zero, this integral diverges: The explicit formula for (i^Q given in appendix 
shows that this constant contains a term that is inversely proportional to 
Mjr. The formula may be viewed as a low energy theorem: The leading term 
in the expansion of (i^Q in powers of the quark masses is determined by 
and F^. At the same time, however, the above discussion shows that this 
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expansion converges only very slowly, because the terms of higher order pick 
up important contributions from the same one loop integrals that determine 
the leading coefficient. 

We emphasize that the infrared singularities seen here are perfectly well 
described by the relativistic loop integrals occurring in our framework. They 
generate problems only if we attempt to replace the loop integrals by their 
chiral expansion, as it is done ab initio in Heavy Baryon Chiral Perturbation 
Theory. The dotted curve in fig.|10| indicates that a decent representation of 
this type may be obtained if (a) the kinematic variables are carefully chosen 
and (b) the polynomial part of the amplitude is dealt with properly. We 
do not elaborate further on the chiral expansion of the loop integrals, but 
now turn to the contributions from the terms of 0{q^) or higher, which our 
representation does not account for. 

19 Higher orders 

The representation of the scattering amplitude to O(g^) incorporates the 
poles and cuts generated by the exchange of one or two stable particles 
(nucleons or pions), but accounts for all other singularities only summarily, 
through their contributions to the effective coupling constants. 

The first resonance in the s- and w-channels is the A(1232). In the t- 
channel, the p(770) plays an analogous role. These states generate poles 
on unphysical sheets of the scattering amplitude. Since we did not include 
corresponding dynamical degrees of freedom in our effective Lagrangian, we 
are in effect replacing these singularities by the first few terms of the Taylor 
series expansion in powers of u and t. 

The p- meson pole occurs at t ~ 30M^. In the region |t| <, lOM^, the 
contribution from this state is well represented by a polynomial, so that the 
description in terms of effective coupling constants is adequate. The range 
of validity of this approximation is about the same as the one encountered 
in TTTT scattering. 

The situation is different for the A-resonance. The curvature caused by 
this state in the low energy region may be described in terms of a low order 
polynomial only crudely. A better description of the effects generated by the 
singularities at s ~ m^, u ~ is obtained from an effective Lagrangian 
that includes the degrees of freedom of the A among the dynamical variables, 
a procedure followed in many papers on baryon chiral perturbation theory. 
For details concerning this approach, we refer to the literature |^S[. In the 
following, we make use of the explicit representation for the tree graph con- 
tribution to the scattering amplitude associated with A exchange given in 
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appendix D of ref. 0. We denote this contribution by D^(z/, t), i?J(z/, t). It 
involves two parameters: the mass and the couphng constant g^. The 
A-resonance starts contributing at O(g^). In a calculation to order the 

(2) 

coupling constants of C}^ enter the loop diagrams, so that part of the A 
contribution is absorbed by the loop graphs with one insertion of a vertex 
from C^^\ 

For the amplitude at t = the situation is simple, because the contri- 
bution of the -couplings to the dispersion integral Df{s) is of 0{q^) and 
numerically negligible. We can thus estimate the remaining A contribution 
by subtracting from the A tree graph the contribution to the polynomial 
part Dp. At threshold, the remainder is given by 

^+ ^ 64 gj ml {{m^ + m^f - Ml) ^^f 



It is formally of but numerically enhanced by the small denominator. 

Inserting the value = 13.0 GeV~^ |28|, we obtain F^D^ = 2 MeV, indicat- 
ing that for a determination of the S-term, corrections of this type are not 
entirely negligible. The corresponding contribution to the scattering length 
is small: Oq^^ = 0.0022 M~^. The effective range of the S-waves or the scat- 
tering lengths of the P- and D-waves involve derivatives of the amplitude 
at threshold, so that the curvature due to the A generates more significant 
effects. 

Note that the singularities generated by the A do not affect the Goldberger- 
Treiman relation. In the case of the low energy theorem that relates the scat- 
tering amplitude at the Cheng-Dashen point to the cx-term, the higher order 
terms due to the A amount to a correction of order M^/(Ma — m,^)^. The 
effect is negligibly small: It increases the value of Aqd by 0.08 MeV. Also, if 
the effective theory is used only in the threshold region, for instance to study 
the properties of pionic atoms the contributions generated by the A are 
adequately described by a polynomial - it suffices to replace the Taylor series 
of the crossing even quantities D^, B^/u, D^/u, in powers of v by an 
expansion in powers of z/^ — M^. The effective value of d^^ relevant in that 
context is the value of at threshold rather than the one at = 0, and 
analogously for the other coefficients occurring in the polynomial part of our 
representation. The straightforward identification of these terms with the 
coefficients of the subthreshold expansion yields an accurate representation 
in the vicinity of z/ = t = 0. It is adequate also at the Cheng-Dashen point, 
but at threshold, the higher order terms of the subthreshold expansion do 
become important, quite irrespective of the contributions generated by the 
branch point that occurs there. 
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20 Values of the coupling constants of C^^^ 



In the representation of the amphtude specified in section |T5|, the couphng 
constants Ci, C2, C3, C4 only enter through the imaginary parts of the disper- 
sion integrals. In section |12|, we have fixed these with the leading coefficients 
of the subthreshold expansion. Numerically, using the values of the sub- 
threshold coefficients of ref. [Q, this leads to 



Ci = -0.6 , C2 = 1.6 , C3 = -3.4 , C4 = 2.0 . (20.1) 

If we use these values to evaluate the scattering amplitude at threshold, 
the result differs significantly from what is observed. For the loop integrals 
to describe the strength of the s- and w-channel threshold singularities, it 
might be better to fix the couplings of there, using the values of the 
amplitude and their first derivatives at threshold as input rather than those 
at z/ = 0. The difference is typical of the higher order effects that our 
calculation neglects and we now briefly discuss it numerically. 

At tree level, the coupling constants ci, . . . , C4 only occur in and . 
The values and flrst derivatives of these amplitudes at threshold are related to 
the scattering lengths of the 5- waves, aQ^_, to their effective ranges, 60+! 
to the scattering lengths of the P- waves, a^_j_ (see appendix 0). Expanding 

the tree graph contributions from C^^^ and C^^^ around threshold, we thus 
obtain a representation of the threshold parameters in terms of the effective 
coupling constants, valid at tree level. Inverting the relation, we arrive at a 
corresponding representation for the coupling constants: 

ci = 7rF2|-(4 + 2a + a2)^ + 6++-3aa++| H-a^Cs , 

2 77^2 

1 + a 

2 7rF2 



c. ^ T^|^ + ^^. + K. + atK8c., (20.2) 



C3 



1 + a 



[{I -a) + a 6+ - (2 + 3 a + 3 a^) a+ - a+ 1 + 16 



C4 = --^ + 47rF2 I - Gr+ + a^-\- (4 - a") , 
4 mjv I 4 mf, ' 



g2 p2 



4 (4 — q;2)2 m: 



With Koch's values for the threshold parameters 0], the above representa- 
tion yields 

Ci = -0.9m-\ C2 = 2.5m-\ C3 = -4.2m-i , 04 = 2.3 m'^ (20.3) 
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Although, algebraically, the formulae ( |12.1| ) and ( p0.2| ) differ only through 
terms of 0{q^), the comparison of the corresponding numerical values in 
( |20.1| ) and ( |20.3D , shows that the difference 



5ci = Ci\ - Ci\ (20.4) 

is quite significant. It mainly arises from the curvature due to the A. In 
leading order of the small scale expansion [^, the shift in the values of C2, 
C3 and C4 generated by this state is given by 



9 (m^ — 



■N J 



\3 



and thus involves the third power of the small energy denominator —m^. 
Inserting the value = 13GeV^^ we obtain 7 ~ 0.23 m~^, so that 



C2 is increased by about one unit, while jcal and C4 pick up one half and 
one quarter of a unit, respectively. In ci, the singularity does not show up 
equally strongly: At leading order in the small scale expansion, we have 
cf = —{m^ — mpf)/mpf'~f ^ — 0.07m~^. The comparison with the numbers 
given above confirms the claim that the A is responsible for the bulk of the 
difference. 

Concerning the t-channel imaginary parts, the role of the coupling con- 
stants ci,C2,C3,C4 and the contributions from the A are discussed in detail 

in ref. 0. These are not sensitive to the manner in which the coupling 

(2) 

constants of £)v are fixed. As an illustration, we consider the quantity 
Act = cr(2M^) — cr(0), relevant for a measurement of the cx-term on the basis 
of TiN data. In ref. ||^, we pinned the coupling constants down with the 
subthreshold expansion and obtained Ao- = 14.0 MeV. If we use the same 
formulae, but fix the values of the coupling constants with the scattering 
lengths, this number changes to Ao- = 15.9 MeV. Both of these values agree 



quite well with the result of the dispersive analysis, A^- = 15.2 MeV [23 . 



21 5'— wave scattering lengths 

At leading order of the chiral perturbation series, the scattering amplitude is 
described by the tree graphs from This leads to Weinberg's prediction 
for the two scattering lengths: 

Numerically, Weinberg's formula yields 4 7r (1 + a) a^^ = 1.14 M~^, remark- 
ably close to the experimental value. 
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The one loop representation of the scattering amphtude accounts for the 
corrections up to and including O^q"^). According to eq. (|15.1|) , the scattering 
lengths are given by 



An{l + a) 
47r(l + a) 



Dt, + 4o + M! rf+ + Mt dtn + D 



pv 



'tt "10 



where and collect the contributions from the dispersion integrals at 
1/ = M^, t = 0. 

In view of the fact that the value of the coupling constant Qt^n plays a sig- 
nificant role, it is not a simple matter to combine the information obtained by 
different authors ||32| . An account of recent work on the low energy parame- 
ters may be found in ref. For the purpose of the present discussion, the 
precise values of the various experimental quantities do not play a significant 
role, however - we simply stick to the KA84 solution |4^. The experimental 
values of the various terms occurring in the above decomposition are listed 
in the table (all quantities in units of M^). 





-Dpi) 




diQ 


d20 




total 


47r (1 -f a) aj^ 

47r (1 + a) a^^ 


-0.15 
0.01 


-1.49 
1.51 


1.17 
-0.17 


0.2 


0.15 
-0.035 


-0.12 
1.32 



The values for Df quoted there represent the full dispersion integral rather 
than the one loop approximation to it, so that the sum of the terms is the 
value of the physical amplitude at threshold. Evaluating these quantities 
in one loop approximation, we instead obtain = 0.08 and D~ = 
— 0.035 M~^: While the value for D~ agrees with the experimental result, 
the one for is too small by a factor of two. Although we are discussing 
small effects here, they do matter at the level of accuracy needed to extract 
the E-term: In the case of Df, the experimental value is = 9.3 MeV, 

while the loop graphs yield = 4.9 MeV. 

We conclude that the one loop representation of chiral perturbation the- 
ory does not cover a sufficiently large kinematic region to serve as a bridge 
between the experimentally accessible region and the Cheng-Dashen point. 
A meaningful extrapolation of the data into the subthreshold region can be 
achieved only by means of dispersive methods - we will outline one such 
method in section First, however, we wish to identify the origin of the 
problem, by studying the difference between the full imaginary parts and 
the one loop representation thereof. We focus on the quantities Df, which 
represent dispersion integrals over total cross sections. 
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22 Total cross section 



The optical theorem relates the imaginary part of the amplitude D{s,t) at 
t = to the total cross section: 

lm,D{s,0) = katot . (22.1) 

To the order of the low energy expansion we are considering here, the scat- 
tering is elastic, so that atot is the integral over the differential cross section, 
given by the square of the scattering amplitude. In the c. m. system, the 
explicit expression may be written in the simple form 

da ml r „ t 



"AT 



4m^ 16 J 

If the chiral perturbation series of the scattering amplitude is truncated 
at 0{q^), the imaginary part of the loop integrals is given by the square of 
the tree graphs generated by (current algebra approximation). Although 
current algebra does describe the 5'-wave scattering lengths remarkably well 
it does not yield a realistic picture for the imaginary parts for two reasons: (a) 
Since we are now dealing with the square of the amplitude, the deficiencies 
of the current algebra approximation become more visible and (b) the effects 
due to the A are more pronounced than in the subthreshold region. This is 
illustrated in fig. where the dotted lines correspond to the current algebra 
approximation. The full lines show the behaviour of the experimental total 
cross section, reconstructed from the results of Koch for the amplitudes 
in the threshold region. 

The discrepancy seen near threshold is directly related to the difference 
between the experimental values of the scattering lengths and the current 
algebra result: There, the cross section is the square of the scattering length, 

a^tot = 47r {a^^ + 2ao_^)^ , a^ot = 4 7r (a^+ - Gq^)^ . 

These expressions are dominated by - the term aj^ vanishes at leading 
order. Although Weinberg's prediction represents a very decent approxima- 
tion, the difference to the observed scattering lengths generates an effect of 
order 25 % in the cross section - this explains why, in the vicinity of threshold, 
the dotted lines are on the low side. 

(2) 

Since the corrections from do not enter the large scattering length 
Og^., the sum of the contributions from and (dot-dashed line) is close 
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Figure 11: Total cross sections. Dotted: Current algebra. Dot-dashed: Tree 



graphs from C^^^ + C^"^ . 



Dashed: CHPT 0(g), all couplings fixed at u = t 



0. Thin: CHPT 0(g ), all couplings fixed at threshold. Thick: Experimental 
values |1T7|. 



to the current algebra approximation at threshold. At higher energies the 
second order terms are however sizeable: The couplings of pick up large 
values because of the presence of the A resonance. Although the corrections 
improve the energy dependence, the one loop representation underestimates 
the cross section, in particular, for I = ^. In the dispersion integral for 
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1 3 

, which involves the difference (72^^^ — cra^oj, the deficiencies happen to 
cancel, so that the one loop approximation for this quantity agrees with the 

1 3 

full result, but for Df, where the imaginary part is given hy tot + '^cr'^ tot, 
the miracle does not happen - this is why the one loop approximation is too 
small by a factor of two in that case. 

The remaining curves show the cross section obtained with the one loop 
representation of the scattering amplitude: The dashed lines result if all 
combinations of the coupling constants Ci,di,ei that enter the amplitudes 
are determined with the subthreshold coefficients, while for the thin lines, 
these are evaluated from the scattering lengths. Since in the second case the 
input parameters are all fixed at threshold, this option yields a somewhat 
better description at higher energies, but in either case the one-loop results 
start deviating from the experimental values already at rather low energies. 
The fact that the dispersive part of the one loop representation is deter- 
mined by the square of the tree level amplitudes rather than by the square of 
the one-loop approximation means that this representation obeys unitarity 
only up to contributions of order q^. The difference between the dash-dotted 
and thin lines shows that the violation of unitarity is numerically quite im- 
portant, particularly in the case of / = |. The figure also shows that the 
perturbative expansion of the amplitude goes out of control immediately 
above threshold: The higher order "corrections" exceed the "leading" order 
contribution. 

The approach used in refs. |^, Q is different. There, the main goal is 
an analysis of the nN scattering amplitude in the region above threshold. 
The coupling constants are determined by fitting unitarized one loop partial 
waves to the data and a rather decent description is obtained - we expect 
that the corresponding results for the total cross section would come quite 
close to the experimental curves in the above figures. The backside of the 
coin is that the representation then unavoidably fails in the subthreshold 
region, so that it is not possible to establish contact with the low energy 
theorems of chiral symmetry. 

23 Analog of the Roy equations 

The preceding discussion makes it evident that the chiral representation of 
the irN scattering amplitude to 0(g'^) can provide a decent approximation 
only in the subthreshold region - a reliable determination of the cr-term from 
the data in the physical region cannot be performed on this basis. Even 
in the case of tttt scattering, where the chiral representation does converge 
rapidly enough for the one-loop approximation to yield a decent description 
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of the amplitude in part of the physical region, the Roy equations provide 
a more accurate framework, which moreover also covers significantly higher 
energies |]TB|, |T^]. In the present section, we briefly outline the steps required 



to extend the Roy equation analysis to the case of nN scattering. 

As shown by Roy [|l^, analyticity and crossing symmetry imply a set of 
integral equations for the partial waves of vrvr scattering. The case of vrvr 
scattering is special in that the s-, t- and w-channels correspond to the same 
physical situation, so that the partial wave decomposition is the same in all 
three channels. The Roy equations represent the real parts of the partial 
waves in terms of integrals over their imaginary parts and two subtraction 
constants. In the elastic region, unitarity then converts this system into a set 
of coupled integral equations. At low energies, the angular momentum barrier 
suppresses the higher partial waves, so that the amplitude is dominated by 
the S- and P- waves. As the pions are spinless and carry isospin 1, Bose 
statistics implies that there are three such waves, to(s), ti{s), tl{s) (the 
upper index denotes the total isospin of the vrvr state, while the lower one 
specifies its total angular momentum). The effects due to the higher partial 
waves can be accounted for in the so-called driving terms. In the elastic 
region, however, these are very small and can just as well be dropped. For a 
detailed discussion of the properties of the Roy equations, we refer to 



Two complications arise when extending this framework to vrA^ scattering: 
The proton carries spin | and, more importantly, there are now two different 
categories of unitarity cuts that matter at low energies: ttN states, relevant 
for the discontinuities in the s- and u-channels, and states of the type vrvr 
or NN, which generate the discontinuities in the t-channel. Accordingly, 
we need to invoke two different partial wave decompositions. We use the 



notation of |28| and denote the s-channel partial waves by f£±{s). For i = 



and i = I, there are altogether 6 such waves: 

fHs) ^ /o^+(^) , fHs) ^ ftis) ftis) = fUs) . (23.3) 

The upper index refers to isospin - it is convenient to use the t-channel isospin 
basis, where /"*" and /~ correspond to J* = and /* = 1, respectively. 

In the t-channel, NN states of angular momentum i carry parity (—1)^"^^. 
They can thus only couple to vrvr states of total angular momentum J = 
In the notation of ref. |2^, the corresponding S- and P- waves are denoted by 
f±{t), where the upper index specifies the total angular momentum J, while 
the lower one labels the two independent partial waves occurring for J > 1. 
In the t-channel, the isospin quantum number is fixed by the total angular 
momentum: J even implies I = 0, J odd corresponds to / = 1. Ignoring 
terms with J > 2, the partial wave decomposition in the isospin even channel 
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contains a single wave, while the isospin odd channel contains two: 



ftit) ^ flit) , /,-(t) ^ fUt) , Uit) ^ flit) . (23.4) 

Altogether, the discontinuities generated by the S- and P-waves thus 
involve 9 functions of a single variable, to be compared with the 3 functions 
to(s), t\{s), tl{s) that occur in the case of nn scattering. A closed system of 
integral equations for these functions is proposed in appendix |K[ It is of the 
same form as the Roy equations: 

^ poo 

I. 1 J Sh 



k=l 
5 



/oo 
dyKr^ix,y)lmf-iy), (23.5) 
- i- 



Sk 



k=l 

(m^ + M,)2 , A; = 1, 2, 3 

4M2 , A; = 4, 5 



The first term on the right represents the partial wave projection of the 
pseudovector Born term, while the second contains the subtraction constants. 
The kernels K^{x^y) are explicitly known kinematic functions. The main 
term on the diagonal is the familiar Cauchy kernel, 

Kfki^, y) = —, — — — + • • ■ 

7r(?/ — X — le) 

In principle, only one subtraction constant is required, but additional 
subtractions may be introduced in order to arrive at a system of equations 
that is less sensitive to the contributions from higher partial waves or higher 
energies. In the case of vrvr scattering, it is advantageous to work with two 
subtraction constants, although, in principle, one of these is superfluous: The 
Olsson sum rule relates a combination of the two subtraction constants to 
an integral over total cross sections. The specific form of the equations pro- 
posed in appendix |^ for the case of t^N scattering involves four subtraction 
constants. One of these can be represented as an integral over total cross 
sections, while the remaining three are to be determined with the experi- 
mental information available at low energies. In particular, the beautiful 
data obtained from pionic atoms |^ subject these constants to a stringent 



constraint. 

The above system of equations may be viewed as a unitarization of the low 
energy representation provided by chiral perturbation theory: The dispersion 
integrals describe the low energy singularities that necessarily accompany the 
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tree graph contributions, on account of unitarity. There are two differences 
to the one loop representation of the unitarity corrections in eqs. ( |16.2|) , 
(|16.3|): (i) That representation accounts for the imaginary parts only to 



first nonleading order of the chiral expansion, while the one underlying the 
above set of integral equations retains the full imaginary parts of the 5- 
and P-waves. (ii) The chiral representation involves an oversubtracted set 
of dispersion integrals. The extra subtractions are needed, because in that 
representation, the energy dependence is expanded in a Taylor series. The 
problem does not occur in the above framework, because unitarity implies 
that the partial waves tend to zero at high energies. 

In the case of tttt scattering, the integral equations can be written in exact 
form, where all of the partial waves are accounted for. The equations for the 
S- and P-waves then contain corrections from the higher angular momenta, 
collected in the driving terms mentioned above. We do not know of a corre- 
sponding set of exact equations for the case of nN scattering - driving terms 
would have to be added to convert (|23.5| ) into an exact system. The size of 



these corrections depends on the number of subtractions made and on the 
range of energies considered. With four subtractions, they should be negligi- 
bly small in the elastic region, because the bulk of the contributions from the 
partial waves with £ >2 is then absorbed in the subtraction constants. The 
available detailed phenomenological information about the imaginary parts 
of the higher waves ^| allows an explicit evaluation of their contributions. 



so that this can be checked. 

Alternative unitarizations are proposed in the literature ||4^ , The 
advantage of the present proposal is that the resulting representation of the 
scattering amplitude does not introduce any fictitious singularities. In partic- 
ular, once the system of equations is solved and the subtraction constants are 
determined, the representation yields a formula for the value of the scattering 
amplitude at the Cheng-Dashen point: 

^ dt, + 2M^rf+ + UMt I dt 



4M2 



t2 {t - 2M2) (4m2 - t) 



oo 



The approach outlined above is closely related to the one used by R. Koch 
T7|, who relies on fixed z/ dispersion relations to study the t-dependence. 



Also, a similar method was used in ref. |^ to extract the value of the cr- 
term from the data available at the time. For more recent work in this 
direction, we refer to , where hyperbolic dispersion relations are used 



to supplement the partial wave relations [E^, for Kn scattering. 
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24 Summary and conclusion 



1. We have derived a representation of the nN scattering amphtude that 
is vahd to the fourth order of the chiral expansion. The infrared singularities 
occurring in that representation are stronger than for tttt scattering. The 
difference also manifests itself in the expansion of the nucleon mass or of the 
cr-term in powers of the quark masses, which contain contributions of O(M^) 
that are quite important numerically. The Goldbcrgcr-Treiman relation and 
the low energy theorem that relates the cx-term to the value of the scattering 
amplitude at the Cheng-Dashen point, on the other hand, do not contain 
infrared singularities, at the order considered. 

2. We have shown that the Mandelstam double spectral function only 
contributes beyond the order of the low energy expansion examined in the 
present paper. The dependence of the box graph on the momentum transfer, 
for instance, can be expanded in powers of t/Am"^ - the cut due to NN 
intermediate states in the t-channel does not manifest itself at finite orders 
of the low energy expansion. Up to and including 0{q'^), the scattering 
amplitude can be written in terms of functions of a single variable, either 
s or t or u. These functions describe the singularities associated with the 
unitarity cuts. In the framework of the effective theory, they arise from the 
one loop graphs. We have given an explicit, dispersive representation for 
these functions, which accounts for all terms to O(g^) and sums up those 
higher order contributions that are generated by the threshold singularities. 

3. The chiral expansion of the dispersion integrals reduces these to al- 
gebraic expressions involving elementary functions such as arccos(— a;/M^). 
On general grounds, the result of this expansion must agree with what is 
obtained in HBCHPT, where the loop integrals are expanded ab initio. That 
expansion, however, is a subtle matter, because the integrals contain infrared 
singularities. In particular, the choice of the kinematic variable to be kept 
fixed plays an important role. While the expansion at fixed u/M^^ does not 
lead to a decent description of the full loop integrals, the one at fixed ojqjMT^ 
does yield an approximate representation for the relativistic integrals, even 
above threshold. 

4. While the expansion at fixed tOq/MT^ solves one problem, it generates 
another. To absorb the divergences of the loop integrals in the coupling 
constants of the effective Lagrangian, the integrals must be expanded in 
powers of M^^ at fixed cj/Af^, so that an expansion of the subtraction terms 
is required, which converges only slowly. Algebraically, the accuracy of the 
representation does not depend on the splitting between the polynomial part 
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and the cut contribution, but numerically, the amplitude is very sensitive to 
this choice. We find it more convenient not to invoke an expansion of the 
integrals - the problem docs not occur in the relativistic representation that 
forms the basis of our work. 

5. Wc find that the dependence of the amplitude on the momentum trans- 
fer t is adequately described by the contributions from the one loop graphs, so 
that the extrapolation from the Cheng-Dashen point to t = does not pose a 
significant problem. The representation obtained for the dependence on the 
variables s and u, however, has a very limited range of validity. Higher order 
effects, such as those due to the A do matter. The one loop representation 
does not provide the bridge needed to connect the value of the amplitude at 
the Cheng-Dashen point to the physical region. 

6. We conclude that dispersive methods are required to obtain a reliable 
description of the scattering amplitude at low energies. With this in mind, 
we propose a system of integral equations that is analogous to the Roy equa- 
tions for TTTT scattering and interrelates the lowest partial wave amplitudes 
associated with the s-, t- and it-channels. The structure of the amplitude 
that underlies this system is very similar to the dispersive representation 
mentioned above, but the constraints imposed on the partial waves by uni- 
tarity are now strictly obeyed. It remains to be seen whether or not this set 
of equations provides a useful framework for the analysis of the low energy 
structure, in particular for a reliable determination of the cr-term. 
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A Notation 



We use the notation of Hohler's "Collection of Pion-Nucleon Scattering For- 
mulas" . (Appendix of 



Kinematics 



P,g 
P ,g 



four momenta of the incoming nucleon and pion 
four momenta of the outgoing nucleon and pion 



s = {p + q)\ t={q- q'Y , u = {P - q'Y , s + t + u = 2M^ + 2m 



4 



E 



s-wl + Ml 



s + ml- Ml 
2^s 



u^l - Ml 
cos 61 = 1 + 

2q^ 



c. m. pion energy 



cm. nucleon energy 



cm. momentum 



scattering angle (c. m. system) 



s-ml- Ml 



2mjv 



k = 



Ml 



lab. energy of the incoming pion 
lab. momentum of the incoming pion 



s — u 
1 

4mA 



t 



UJ 



{t - 2Ml 



4:171^ 

2\ 



Isospin 

We denote the amplitudes for the reactions Tr^p n^p by A±, the charge 
exchange reaction amplitude ir^p Tc^n by Aq. The amplitudes of the 
isospin eigenstates are denoted by A^, J = |, |. The relations to the isospin 
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odd and even amplitudes A"^ are 

A+^Ai = A+-A-, A_ = 1{2A^ + A^) = A+ + A- 

Ao = ^(Ai - A5) = -V2A- , = i(3A_ - A+) ^A+ + 2A- 

A+ ^ 1{A^ + 2 aI) , A-^l{A-2-A^ 

Partial waves 

For the partial wave decomposition, the amplitude 

is more convenient to work with than A ox D. 
In the s-channel, the decomposition reads 



C = '^^Y.^^+l)h+{s){{E + m^)P^{z)-{E-m^)P^+^{z)] 
-'^^Y.lf,_{s) {{E + m,)P,{z) -{E- m^)Pt-,{z)} , 



4 



B = ^5^^+(s){-(^ + m.)P;(;.) + (E-m.)P;+i(^)} (A.l) 



1=1 

oo 



The i-channel partial waves are defined by 

J=0 

oo 



^ = , (A.2) 

j=i V ■^(■^ + ^) 
Zp-q- 
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Threshold expansion 

We write the threshold expansion of the amphtude X G {-D"*", D~ , B~^, B~} 
in the form 

RcX{q^ t) = Xoo + Xio + X^.t + X^o + Xn qH + Xo2t^ + ... 

The coefficients X„„, of this expansion are related to the scattering lengths 
and effective ranges. For the threshold expansion of the real part of the 
partial wave amplitudes 

Ti^ = qfi± = ^{v±e'''^^-1) 
zi 

we use the notation 

ReT;± = g''+^(a,± + + q± + . . . ) ■ 

The lowest coefficients of the threshold expansion of the amphtudes D and 
B are obtained from the expansion of the partial wave amplitudes via 

Doo = 47r (1 + a) ao+ 

Dio =47r(l + a) |-^^ + 6o+ + ai- + 2ai+j 
L>oi =2 7r|^ + ai_ + ai+ (2 + 3q;)| 

^ 2 71 (1 + a) ( ,^ 2\ "^0+ , r. 2 

D20 = 5 — i - (1 - a + a ) ^ „ „ + 60+ + 2 q; co+ + ai_ 

+ 2 a 6i_ + 2 ai_|_ + 4 a bi+ + 4 a 02- + 60; a2+ 1 

+ 4(2 + 3a)m2 61+ + 12 (2 + a) 02- + 12 {3 + Aa)mla2+^ 

D02 = 2— V {'^1+ + 4m^a2- + 2 (3 + 5 a) 02+ 1 

5oo = 8 TT m^v I + ai_ - ai+ 1 

2 TT f ao+ 9 

-Bio = — i - — 5- + O0+ + ai- + 4 m^^ 6i_ + 2 ai+ 
m^r l 4mj^ 

- 4m^6i+ + 12m^a2- - 12m^a2+| 

3 TT f 2 i 2 1 

-Doi = — i ai+ + 4 a2_ - 4 m^^ 02+ ^ . 
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B Off-shell TrAT-amplitude 



As discussed in the text, the off-shell amplitude of the effective theory is 
without physical meaning. In the underlying theory, however, an unambigu- 
ous off-shell extrapolation can be constructed from correlators of currents 
with the appropriate quantum numbers. For the pion field, the operator of 
lowest dimension is the pseudoscalar density P"(x) 



To generate nucleons, we need a three-quark operator with the correct quan- 
tum numbers. At lowest dimension, there are two independent such opera- 
tors: 



In qa{x) = {ua{x) , da{x)) we have collected the quarks with color index a. 
The brackets indicate a sum over Dirac and flavor indices, ir^ stands for the 
antisymmetric tensor in SU(2) flavor space and C = —i^oj2 is the charge 
conjugation matrix. These operators couple to the proton with coupling 
strength Gi,2- 



Since these operators carry anomalous dimension, their matrix elements de- 
pend on the running scale of QCD, but arc otherwise free from ambiguities. 
Their Green functions contain contact terms that make it difficult to han- 
dle the corresponding generating functional - the perturbation generated by 
the source terms is not renormalizable. We may, however, restrict ourselves 
to Green functions involving only two nucleon fields, so that a finite num- 
ber of counter terms suffices. The pion nucleon scattering amplitude can be 
extracted e. g. from the Green function 



Under a chiral transformation q{x) — >• -|- 75) Vr -|- |(1 — 'y5)VL}q{x) the 
operator Ni transforms as 



P"(x) = q{x) i^^T" q{x) , 1,2,3 
which couples to the pion with strength G,,, 

(0| P'^(O) I TT^) = 6-^(5"^ . 




{0\N,,2{x)\P,s)^GiMu{P,s) 



(0|T Ni(x')P"'(?/')^"(?/)Ni(x) |0). 



Ni ^ {^(1 + 75) ^ii - ^(1 - 75)^l}Ni 



63 



In our effective theory, it is easy to identify an object tliat transforms in 
tfie same way, namely \l/ = uipR — ipi. If we would add a source term 
f^i Ni + Ni T] to the Lagrangian of QCD, this source would couple to the field 
in our effective Lagrangian and we would be able to compute the above 
Green function with an extended version of our effective Lagrangian. Since 
the two fields '^{x) and ip{x) do generate the same on-shell amplitudes we 
can also use the field ip{x) to extract these, as we did in our calculation. 



C One-loop integrals 

In this appendix, we list the integrals that occur in the evaluation of the 
scattering amplitude to one loop. Throughout, we use infrared regulariza- 
tion (for a detailed discussion of that method, we refer to 0). We perform 
the tensor decomposition of the integrals with respect to the sums and the 
differences of the external momenta, 

= (P + qY = {P' + q'Y 
QA* = {P' + PY 

= {q - qY = - P'Y ■ 

and put all external legs on their mass shell, 

P2 = p'2 ^ ^2 ^ ^2 ^ ^/2 ^ ^2 _ 

The coefficients represent Lorentz invariant functions of the Mandelstam vari- 
ables. The squares of the vectors S, A and Q are given by 

= s , A^ = t , = Am'^-t 



1 meson: A^r = Jio 



1 nucleon: A 



N — J-01 



I f dfik 1 ^ 



" i J J {2tiY - 

Note that integrals formed exclusively with nucleon propagators vanish in 
infrared regularization: /o„ = 0. 
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2 mesons: J = I20 



{J, J^ jn = - / 
^ J I 



{2ttY (M2 - A;2)(M2 - {k - A)^) 
1 meson, 1 nucleon: / = In 

d'^k {1 , kf" , kf" k" } 



^ J I 

2 mesons, 1 nucleon: 



(27r)'^ (M2-A;2)(m2-(S-A;)^ 



{-^21 , -^21" 5 hi"} — - 



1 f d'^k {1, A;'', fc'^fc'^} 



z J I {2ttY (M2 - P) (M2 - (A: - A)2) (m2 - (P - k)^) 

I,r = 9'''lil\t) + Q'-Q'^lif + A'^A^J^f (t) + (Af^Q^ + Q^^An'^lS\t) , 
1 meson, 2 nucleons: /12, /a, Ib- 

d'^k {1 , kf" , k^' k" } 



^ J I 



i J I {2tiY (M2 - (m2 - (Pi - kf) (m2 - (P2 - kf) 
With Q = Pi + P2 and A = Pi — P2, the tensorial decomposition reads 

i^^'' = gt'^if^ + g'^g'^/fj + A^A^jfJ + (A^g^ + g^A^)42^ . 

Several different topologies give rise to this integral: the graphs (c), (d), (g), 
(h) and (m) in fig. ^. The values of Pi and P2 differ from graph to graph. 
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In all cases, however, one of the two momenta is on the mass shell - we set 
Pf = m^. The coefficients of the tensor decomposition then only depend on 
two variables: 

/S) = jS)(.,t), P! = m\ Pi = s, t = {P,-P,f. 

In the case of the graph (m), we have Pi = P, P2 = P', so that we are dealing 
with the special case s = w? . In the case of (c), (d), (g) and (h) on the other 
hand, the difference Pi — P2 represents the momentum of an on-shell pion, so 
that t = M^. To simplify the notation, we introduce invariant functions that 
correspond to these two special cases and only depend on a single variable: 

lA{t) = Ii2{m\t) , lBis) = Iu{s,M^) . 

The coefficients of the decomposition of the corresponding tensorial integrals 
are 

I Tit) = 9^''lf{t) + Q^Q^'lfit) + A^A^jjf)(t) 

1 meson, 3 nucleons: 

1 r d'^k {1, k^'} 

~i Ji{2TrY (M2 - A;2)(m2 - (P - ky){m^ - (S - k f){m^ - {P' - k f) 
I^,{s, t) = {s, t) + ( A + 2qr /i? {s, t) 

D Results for the loop graphs 

In the present appendix, we list the contributions from the various topologies 
shown in fig. |^. Note that the amplitude in addition contains the contribu- 
tions obtained from the graphs (a)-(i) and (n)-(s) through crossing, 

At{s,t) ± At{u,t) 

P^(s, t)TBt{u,t) G E {a, . . . ,i,n, . . . , s) 
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D.l Loop graphs of C^^^ 

The contribution of the loops formed exclusively with vertices of the lowest 
order Lagrangian C^^'^ starts at order 0{q^). We give the exact result for the 
C^^^ loops in terms of the amplitudes A{s,t,u), B{s,t,u), since the corre- 
sponding expressions can be written in a more compact form. Note that the 
amplitude A{s,t,u) does not obey the chiral power counting rules. 
We use the following abbreviation: 

F{s) = -2 I{s) + (s - m^) I^^\s) . 

graphs a+b 

4 — 4+ R~ — R+ 

^ah — -^ab ^ab ~ ^ab 

graphs c+d 

= ^{ - 2 A. + (. - m^)I^^\s) - Sm^ (m^ - (. - m^) I^is)) } 
I is) + {s- m') - W /«(.)) 

A, + (3m2 + s) (m^ Ib{s) - {s - m') I^^\s)')] } 



8F4 

4m^ 



+ 



s — m? 



'^cd '^cd '^cd '^cd 



graph e 
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graph f 



s — m 



graphs g+h 



s - m') I - 2 /(s) + /(^) (s) + 8 /^'^ (s) } 



gh 2 ^4 



graph i 



3g\rn 
16 



|2 (/(m^) - I{s)) + {s- m") (2 l{s) + /(^^ (s)) 



+ 8 [ - 7a(^) + 4 7^^^ (i) - (s - -u) 7^^^ {t) 
+ (75(.) - 7«(5) - 7g)(.)) - K + 3.) 4'^^) - - 
+ 32m^(s-m2)7}3^(s,i)| 



5+ = I (3 + s) 7(s) + 4 771^7(1) (m^) - (m^ + s) 7^^) (s) 

+ 4m' [M2 7^(i) -2 7i')(t) + 2M2 7B(s) -2(3m2 + s)7^')(s) 
+ 2 (m^ - s) 7^')(s)] + 16 [ - 7i3(s, t) + 2 (s - m") 7g)(s, t)] } 



3^^ 



b: 



1 o- 
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graph k 



graph 1 

A" = f^{'2M'H'n') - (M'' -2t) - 4m'4\\t)) } B+ = 

AT = - u) /<?(t) BT = 7«(t) + 4mV<?(t)} 

graph m 



A-^ = -'^{s-u)lf{t) 

= - {^'T - 4m^ (/(^)(m^) + U{t) - 2 if (t)) } 

graphs n - V 



w 



A". = 



At = 



^ ^ , 2m 1 



2m 1 



P'- fMm2-s 2m 



Bf^O 



Btu = 



B. 



tu 



Atu 
2m 



S+ = 



b: 



5A^ 
8F4 
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D.2 Loop graphs of C^"^^ 

The loops which, besides vertices from CS^^ ^ also involve one vertex from CS^^ ^ 
start contributing at order . As exemplified in section (^), the fact that we 
only need the leading order of the expansion of these diagrams simplifies the 
calculation considerably. The results are given in terms of the amplitudes 
B{s,t)} to order {O(g^), 0{q^)}. We define 

(m^ - sf + 2m^{t-2 M^) 



graphs a+b 



D 



ab 



Qit) Bt, Bt, = T;^^^f^^(c3 - C4 (d - 2)) /(^)(.) 



{m? - s) 



graph f 



D^ = B^ = B- = Q 



Df = ^{^ciM^ (m^ - s) J(s) - [8ci - (c2 + C3) (m' - sf] I^^^ (s) 

- 2 (C2 + C3) (m=^ - s) + m^I^^\s)) 



graphs g+h 

= -^ll^ {4<=3<3W4"(s) + Q(0) [- 8c, Af' /<?',,) 

-4(C2 + C3)K-S)4^)(^)]} 



+ - 

gh 




+ _ 
gh 


Bgh 






~ — 




gh 




— 


8 C3 g\ 


gh ~ 
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graph k 



Dt - ^{(^' - (4ci + C3 (t - 2M^)) J(t) 



+ 2(4ciM'- -C3 (M' + 2t))A^-2c2t(2i-M2) J«(i)} 



graph m 



= (c2 - uf - 32 ci + 8 C3 m^(2 - i)) 



S+ = = 

m p4 

graphs n+o 



[{d-5)l^^\t)-tl^\t)} 



D~ ^ D'^ B~ = 

no no no no 



graph s 



Dt - Qit) Bt Bt = ^f^, (7 - 2 + ^3) A. 



d; = b: = Bt 



graph V 
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E Renormalization of the effective couplings 



The renormalization of the couphng constants ^3 and I4 that occur in the 
mesonic part of the effective Lagrangian is discussed in pO[. These constants 
enter our amphtude when expressing the bare quantities M, F in terms of 
their physical values M^, F.^ and in the wave function renormalization of the 



pion field. The renormalization of the couplings of is given m 
the HBCHPT formalism is used to extract the divergent part of the one loop 
functional. In infrared regularization, the amplitude is of the form M'^f, 
where / is a function whose chiral expansion only contains integer powers of 
the chiral expansion parameter. Hence the poles at d = 4 always appear in 
the combination {d — 4)~^ + In M//i. The same is true for HBCHPT, because 
the nonrelativistic loop integrals do not involve the nucleon mass scale. The 
coefficients of the chiral logarithms must be the same in the two formulations 
of the effective theory, because these are of physical significance and cannot 
depend on the regularization scheme. We conclude that the divergences 
encountered in our approach are the same as those occurring in HBCHPT. 

As a side remark, we note that in dimensional regularization, the renor- 
malization of the effective coupling constants is different. In that approach, 
the infrared regular part of the amplitude is retained. This part also con- 
tains poles at (i = 4, but these are not accompanied by a chiral logarithm. 
Instead the pole terms occur in the combination {d — 4)"^ + Inm/yU. Also, 
some of the divergences contained in the regular and singular parts cancel - 
the dimensionally regulated amplitude contains chiral logarithms without an 
associated pole at d = 4. 

We define the renormalized coupling constants by 



21 , where 



A 



di 



d-i 



1 



ln47r + r'fl) + 1 



(47r)2 yd -A 
For the pion couplings /s, I4, the coefficients are 
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1 

"2 ' 
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2 . 



Those relevant for rfi, ^2, 



5i 

5i4 



S2 

5l5 



12 T?. 3 



may be taken from ref. 
2 



12 







S3 
5w 



^^ = k + YA9' 

5l8 = 



These expressions differ from those given in [|^ because these authors intro- 
duce additional equation of motion terms in the Lagrangian, to arrive at a 
finite scattering amplitude also off the mass shell. 
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Finally, we list the coefficients occurring in the renormalization of the 
coupling constants Ci, 62, ... that enter the amplitude at order O(g^): 

ei = I - 1 (8 ci - C2 - 4 C3) m 

£2 = -f fl-^ - I (c2 + 6 C3) m 

€3 = — 1 — 3 (7^ — y (yf"^ — 8 Ci m — C2 m + 4 C3 m 

e4 = 10 + 12 / + f / + 8 C2 m 

65 = | + |5'^ + y5''^-8cim + YC2m + 5c3m 

eg = -12-8(72-8(7^ 

£7- -3 -3^ --g 

es = -5'^ ~ I fi'^ ~ I (C2 + 6 C3) m 

eg = 1 + ^ + 4 C4 m + (13 + 24 C4 m ) 
eio = -f ^?'-|/ (3 + 8c4m) 
^11 = I ~" 6 5'^ - 1 + 1 C4 

One readily checks that these renormalizations ensure a finite result for the 
quantities m^r, (yf^^r, M^, F^, for which we listed the explicit expressions in 
section ^. The same holds for the coefficients occurring in the polynomial 
part of our representation (see appendix |T^), so that the scattering amplitude 
approaches a finite limit when the cutoff is removed. 

F Subthreshold coefficients 

In this appendix, we list the explicit results obtained for the expansion of 
the subthreshold coefficients in powers of the quark masses. These relations 
specify the polynomial part of the dispersive representation (|15.1|) in terms of 
the effective coupling constants. We express the result in terms of the physical 
pion mass. To simplify the formulae, we set the renormalization scale /i equal 
to and indicate this with a tilde on the renormalized couplings: 

di = 

ei = 

Note that the chiral 



ci.--^A. = ciI(/i)-^^-^ 







In 








In 









p2 - 16 7r2F2 

logarithms are then contained in the coupling constants. 
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Amplitude D+ 



_ 2{2c,-c,)Ml , gl{?, + %gl)Ml 
F3 ^ 647rFi 



2c2 (4 + 5^1) 



32 7rF,4 



64 16CiC2 1 + ^1 197 



_ _ (77 + 48ffl)M. 
°^ 7687rFi 



TT TT 



. , 65 , 52ci-C2-32c3 (7l(47 + 66(?i: 



F2 192 7r2F4 3847r2F4mjv 



12 + 5(71 , ee , U + 10 gl , 173^j 



2° 192 TT F2 ^ 24 280 tt^ 



11 C/i 174 A/f _ 172 oi:; 174 ™ 



_ 193 gl 29 gj 



A 



15360 TrFiM^ 8Fjm2 4807r2F^m^ Ut:'^ F^m^ 

19 cx 7c2 7c3 



4807r2Fi 6407r2Fi 807r2Fi 



Amplitude D' 



1 A{d^ + d2 + 2d,)Ml g\Ml 

dan = 77-^ H 7^ \- 



00 2F2^ F2 ^ 487r2Fi 



TT TT 



1287rFim„ 47rFi 
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44 _ 15 + 7^1 (168 + 138(71 + 85^1) 

F2 2407r2F4 768 7rF4m^ 



(8ci + 8c2 + 8c3 + 5c3 5l-5c4^i) 



16 7rF4 



^ 2 (Ji + 4) _ l + 7ffi + 2gl (12 + 53gl + 24ffl)M. 
°' 192 7r2F4 ^ 3847rF4m^ 

_ (C3 - C4) ^1 
87rF4 

Amplitude S+ 



, + _ 4m;v (rfi4 - dis) _ g\mi, gl (8 + 7gl) _ {03 - C4) gl tjIj, 



Amplitude B" 



gl (3 + 2gl + 9c4mJ 
12 7r2F^ 



fif^m^v ^ eipm^r ^ t/^ (25 + 36^1 + 80c4mjv) 



^° 32F4M^7r F2 1207r2F4 



_^eiim^ l-Qgl-Agl + Aam^ 
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G t- channel partial waves 



Only the two vrvr partial waves 

Of — f — A M"^ 

327rF2 ^"""^ 967rF2 

enter the t-channel unitarity relation ( |14.5| ) at the one loop level. The imag- 
inary parts of the amplitudes D, B can thus be constructed from the three 
lowest partial waves /°, of the tree level vrA^ amplitude: 

Imt5+ = , ImtE^ = 6 TT Im . 

The Born term contribution to the relevant t-channel partial waves reads 

TT \ (t-2M2)/t2 24 m2 
91n r (1 + «:^) /(«:)- 1 1 



v/2 7r I (t-2M2)/t2 12 m; 

while the one generated by is linear in t: 

2 

flit) = ——{-24M' ci + {4M'-t)c2 + Q (2 - t) Cs}+0{q') , 
24 7ri' ^ ^ 

^ m+_C4t .1 . ^ 1 + 4 C4 

^" 247rF2 ' l2v^7rF2 ' 

The function /(k) stands for 



_ arctauK _ ^/(t - 4M2) (4m2 - t) 

Note that the chiral expansion of this function does not converge at t = AMI: 
Since k counts as a quantity of 0(g~^), the arctan is replaced by its Taylor 
series in inverse powers of k. The point t = AM^ corresponds to k = 
and is thus outside the radius of convergence of that series. The infrared 
singularities occurring in the scalar form factor in the vicinity of the t-channel 
threshold are discussed in detail in ref. - the scattering amplitude exhibits 
essentially the same structure there. 
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H Imaginary parts of the loop integrals 



In this section we give the exact low energy imaginary part of the loop in- 
tegrals defined in Appendix 0. We can distinguish three different types of 
singularities: 

1. Singularities in the low energy region, 

2. Physical singularities in the high energy region, 

3. Unphysical singularities generated by the regularization. 

The singularities of the first kind are those generated by intermediate states 
that involve at most one nucleon. They are identical in dimensional and 
infrared regularization. The singularities of type (2) and (3) are absent at 
any finite order of the chiral expansion. An example for type (2) is the t- 
channel cut due to NN intermediate states in the triangle integral with two 
nucleon- and one pion-propagator. This integral develops an imaginary part 
for t > 4 m^, which can be obtained from the elastic unitarity relation in the t- 
channel. In the chiral expansion, the corresponding contribution is converted 
into a Taylor series in t. In infrared regularization, the loop integrals contain 
further singularities, which are not related to physical intermediate states. 
The infrared part of the self energy integral /(s), for instance, has a cut 
for s < and a pole at s = 0. Singularities of this type do not occur 
in dimensional regularization. In the following, we only list the imaginary 
parts of type (1) - only these are relevant for the low energy structure of the 
scattering amplitude. 

1 meson, 1 nucleon 

The dimensionally regularized self energy integral has a right hand cut start- 
ing at s = (m + My. Its discontinuity is given by 

Im H{s) = e{s - s+) , s± = {m± Mf 

16 TT S 

and p{s) = a/(s — s+) (s — s_) = 2 m \J uj'^ — . 

The regular part of -ff(s) has a pole at s = and a cut for s < 0. The 
imaginary part of R is 

The infrared part is given by the difference /(s) = Hi^s) — R{s) and hence 
has a cut on the left as well as on the right. The left hand cut lies outside the 
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low energy region and is absent to any order of the chiral expansion. Note 
that the right hand cut is the same in both regularizations. The imaginary 
part of the coefficients of the vector and tensor integrals reads 

2 s Im /(^^ (s) = (s - + M^) Im /(s) 
12slm/(2)(s) = -p{syiml{s) 

3s'^lml^^\s) = (^{s - + M'^Y - Im/(s) 
2 mesons 



1 / 4A//"2 



4 Im J(2) (t) = Im J{t) , 12 i Im (t) = - 4 M^) Im J{t) . 
2 mesons, 1 nucleon 

For t < Am^ the imaginary part of the scalar triangle integral is given by 



, , , , e(t- 4M2 J(t-4M^)(4m^ -t) 

Imhiit) = \ ' : arctan-^^^ ^^-^ . 

^ ' 87rVt(4M2-t) t - 2M2 

Im/2i^(t) can be expressed in terms of Im/2i(t) and Im J(t): 

2 g2 Im it) = (2 - t) Im I21 {t) + 2 Im J(t) , 

8 Im (i) = (4 - Im hi (t) + 2 {t-2 M^) Im {t) , 
8Q^ImI^f{t) = 6 (2M2 - lml^^^(t) + {t -AM^) Im/aiW , 
8tlm4f{t) = 2 (2M2 - Im7^i^(t) + (3i - 4M2) Im/siW , 

where = 4m^ — 
1 meson, 2 nucleons 

We restrict ourselves to the two special cases 

lA(t) = /l2(m^ , Ib(s) = 7i2(s, M^) . 



78 



The functions P^'it) 



are analytic in the low energy region, 



Im/7(i) = 



for i < 4 . 



For s > the absorptive part of Ib{s) is given by 



Im I b{s) 




The imaginary part of the coefficients of the tensorial decomposition can be 
obtained from the relations 



2 p{s f Im 1)^' (s) = (m^ - - s) 

X ((m^ - s + 2 M^) Im Ib{s) + Im I{s)) , 

2 p{sf Im (s) = (m^ - s) (3 + s - 3 M^) Im 

+ (m^ + 3 s - M^) Im I{s) , 

4 Im ^ (s) = 2 Im /b (s) - (2 m u; + 3 M^) Im 7^^^ (s) 

-{2mu; + M^) Im/f^(s) , 

4 p(s)2 Im (s) = 4 Im 4^ (s) - (2 m + 2 M^) Im /^^^ (s) 

+ 4 (2m2u;2 + M2mu;-M^) Im4^(s) , 

4 p(s)2 Im 7^^^ (s) = -4 (2 + m a; - M^) {2mu; + M^) Im (s) 
+ 4 (4m^ + 4ma; + M2) Im7^^^(s) + 2 (2m2 + 3mu; + M^) Im7W(s) , 

4 Im 7|^^ (s) =2 {2muj + M^) (^Im 7^^^ (s) - 2 Im 7^^^ (s)) -Im 7^^^ (s) . 
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1 meson, 3 nucleons 

The absorptive part of the box integral Ji3(s,t) is a function of t. For t < 
4m^, it is given by 



lmii3(s, t) = arcsmh 



C(s) = {2mLU-M^) (m^ s - (m^ - M^) 

4 {p{sy + st) Im/J5^(s,t) = (m^-M^-s) |2Im/B(s) 

+ {2m^ -2s + AM^ - t) ImJi3(s,t)| 

4 (p(s)2 + St) lml[l\s,t) = 2 (m^ - + s) Im/B(s) 

+ ((m^-s) (4m2-2M2-t) - (2M2-t)) Im/i3(s,t) 

Note, that the tensor integrals are not singular at p(s)^ + st = 0. This can 
e.g. be seen by using the relation 



^ s 

To the accuracy of our calculation, we need only the first two terms of the 
Taylor expansion of these integrals in t (see section p.4.3|) . 



I Low energy expansion of the loop integrals 

We give the explicit form of the coefficients of the tensor decomposition of 
the loop integrals to next-to-leading order. The result is expressed in terms 
of the dimensionless variables 

^ s - - t , M 

U = , r = — - , and a = — . 

2mM M2 ' m 
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It involves the three functions 



Vl — arccos(— f2) 




The function f{fl) is associated with the scalar self energy integral and g{T) 
is related to the triangle integral with two meson and a single nucleon prop- 
agator. Note that (?(r) contains arbitrarily high orders of the expansion 
parameter a. It is constructed in such a way that the representations for 
I^{t) also cover the region around r = 4, where the chiral expansion of 



these integrals breaks down (see section |T8| and ref 



2 mesons: J = L 



20 



At) = J{r) 



16 7r2 



2 



J«(t) 



r - 4 - 
12 r ^ ' 



1 



576 7r2 



r o 



4 64 7r^ 2 r 



1 meson, 1 nucleon: / = In 

I{s) = -a (l-2an) f(n) + ~ - 2a (n + a-2an') A^ + 0(a3) 

16 vr^ 

J(l) (S) = + « _ 4 « ^ 1 )_ 

16 vr^ 
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3 14:4: TT'' 

- J (3 Q - 2 + 3 a - 12 a + 8 a Q^) + 0(q;^) 



3 3 O 

7(3) (s) = ^ (l-4Q^-12an+24aQ^) (7Q^-3+18q;Q-30q;Q3) 
3 72 TT^ 

+ -a^ {3n-4n^ + 3a-24aQ^ + 24an'^) + 0{a^) 
3 

2 mesons, 1 nucleon 

hi {t) = (r ) + a A^) + O (a) 



iS'it) = -^{a (2 - r) 5(r) + 2j(r) - - 4 A J + 0( 

8m^ L 8 7r^ J 



= -{2 (4 - r) ^(r) + a (r - 2) J(r) + ^ 

+ 4a(4-r) X^^ + 0{a^) 

(0 ^ ^{2 « (r - 4, + 3 (2 - r, J(r) - ^ili^Zj^ 

+ 8Q;2(r-4) A^} + 0(q;^) 

4t'(*) = T^{2 (3 r - 4) + (2 - r) a J(r) + 

A^} +0(q;) 

1 meson, 2 nucleons 



lu {s, t) = - (2 O - a - 2 a Q^) (1 - a A. 

+ 9 ,^j 27rQ + 2Q'-7ra-2aQ + 2QiQ^)+0(a') 
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a 



;(l81]2 + 37ra-21aO-8aO^) 



1152m2 7r2Q 

^a(60 + 9Q;- 16q;Q^) A^ + 0(q;^) 

12 ^ ' 



4? ^) = (n^ - 1) (2 Q - a - 6 a Q^) /(Q) 



12^2 

„2 



+ ^i77rT7|2 0(37r + 6f7-5f)3) -a(37r + 6f7 + 13fl3- 180^)1 

576 TT^ t V /V y J 

(3 -2Q2 -8Q;Q + 6Q;fi^) + O(a^) 

6 



(5, t) = 4^^^ {2 O - 8 03 _ ^ (1 ^ 14 _ 3g I ^ 

■|2n(37r + 6Q- 14Q=^) -Q;(37r + 6Q + 22Q3-72Q^)} 



a2 



2304 m2 Q? 



+ -— r(3-402 - 16aO + 18af^3) + O(a^) 
12 



^^?(^'^)--OT(^ + ^^^)^^^)+1152^(3^ 
1 meson, 3 nucleons 



4 q; 
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ii?is,t) = -:^2-^(2 + n') m) -T7^,K 

+ ^^.o \ 203 (6 + 12 - + 0(a) 
1152 m^TT"^ ir ^ 

J Low energy expansion of the amplitude 



The representation of the amphtude given in sections |T5| and ^ consists of 
the Born term, a set of dispersive contributions that account for the cuts in 
the S-, t- and w-channels and a polynomial. The dispersive part is described 
by the functions Df(s), . . . , Bf{t), which can be expressed as integrals over 
the corresponding imaginary parts, according to ( |16.2|) and ( p.6.3| ). As dis- 



cussed in section |T8|, the first few terms of the chiral expansion of these 
integrals can be evaluated in closed form. In the present appendix, we give 
the corresponding explicit expressions. 

The chiral expansion of the scattering amplitude involves a choice of kine- 



matic variables to be kept fixed. As discussed in section |T8|, the convergence 
of the expansion improves considerably if the lab. energy Q = cu/M^ is re- 
placed by the cm. energy Qg = Uq/M^^. At the same time, however, the 
change of variables generates bookkeeping problems, because polynomials 
in VLq do not represent polynomials in z/, t. For this reason, we work with 
the variables a;,t in the s-channel and use for the t-channel. The cor- 
responding expansion at fixed VLg is readily obtained from the expressions 
given below: It suffices to express uj in terms of VLq and to expand the result 
in powers of a at fixed VLq. 

When combined with the chirally expanded polynomial part specified 
in appendix ^ the expressions given below can be used to work out the 
chiral expansion of the scattering lengths and effective range parameters, 
for instance. As discussed in section the expansion of observables that 
involve derivatives of the amplitude converges only extremely slowly because 
of the strong infrared singularities generated by the s- and w-channel cuts. 
For numerical analysis, we strongly recommend the use of the dispersive 



representation specified in sections O and 16 
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The functions £>f (s), Z)J(s), Bf (s) 

For all of the expanded one-loop graphs that contain a cut in the s-channel, 
the discontinuity across the cut is given by a polynomial in uj and t times 
the factor ^Juj"^ — M'^/uj'^ , where uj is the pion lab momentum. The factor 



^Juj'^ — represents the imaginary part of the elementary function 



arccos 



UJ 



which does develop the required square root discontinuity: Approaching the 
real axis from above, we have 



1 - 



Ml 



-arccosh(— jg-) 
arccos (—-^) 
arccosh(jg-) - 



ITT 



for UJ < -M^ 
for -M^<uj < 
for UJ > Mt, 



The imaginary part determines the real part up to a polynomial. In fact, in 



the normalization specified in section |T6|, the real part is fully determined by 
the imaginary part, through the condition that the Taylor series expansion of 
the functions Df{s), Df{s), 5f (s) in powers of u does not contribute to the 
subthreshold coefficients listed in eq. ( |16.4|) . It therefore suffices to list the 
contributions proportional to f{uj), which we denote by Df{s), Df{s), Bf{s), 
respectively. The explicit expressions read: 



,4 



{At + g\At + g\At}f{uj) 



Uuj'^i-mj^uj + Su^ - M;) 
24 uj^ {uJ^-Ml) 

^{uj^ -Mlf {-m^uj + ?,uj'^ + Ml) 



Dtis) 



Btis) 



{u^ - Ml) 



3 + uj + 2uj^ + Ml) 

- 8 (C3 - Ci)m^uj^}f{uj) 
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Dlis) = 12 .;3 + 9l A, + 9^ A3 + A4 }/(^) 

= 6 tu^(-m^ u + Suj'^ - Ml) 
A2 =6c<;2(2c<;^-2cj2M^ + M^) 
A3 = 2 (-m^ + 3 cu' + M^) (cu^ - M^)^ 
AJ = -48 cu^ { (C2 + C3)cj2 - 2 ci } 

- 16 g\ (C3 - C4) cu^ (cu^ - M^)^ 



^2(g)= i2ir4^^^3 +^^(-m^o; + 3a; + M J 

- 8 (c3 - C4)mjvu;2}/(u;) 
{s) = ^^^^3 -^ {2^^ (-m^c. + 2 a;^ + M,^) + 3 c.^ 

+ 8 C4 mjv t<;^}/(t<;) 

Note that the expressions contain fictitious singularities at = 0. These drop 
out when evaluating the functions Df{s), Df-(s), Bf{s), which are obtained 
by simply removing the leading terms in the expansion in powers of u: 

4 

Dfis) = J2df^^u;- + Dfis), 

n=—3 
2 

n=—3 
2 

Btis) = J2^tn^'' + Btis). 



n=— 3 



The functions Df{t), B^{t) 

Only the diagrams belonging to the topologies (k) and (1) have a branch 
point at t = 4M^. Note that in the case of (k), the two-pion- vertices from 
C^M^ are relevant as well as those from C^^\ which are proportional to ci, 
C2, C3, C4. Explicit representations for the contributions from these graphs 
are given in appendix 0. Those of type (k) may be expressed in terms 
of the familiar loop integral J{t). The low energy structure of the vertex 
diagrams (1) is discussed in detail in 0. To the order we are considering 
here, the corresponding integrals may be expressed in terms of J{t) and of 
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the function g{t) which is the expanded version of the triangle integral l2i{t) 
with two mesons and one nucleon. In the t-channel, we thus encounter two 
functions that play the same role as the quantity f{uj) used for the s-channel. 
Again, it suffices to list the contributions proportional to these functions (the 
explicit expressions for J{t) and g{t) are given in section |l^): 

Dtit) = (24 c. Ml + c,{t-A Ml) 

Bt{t) = 



^3" (t) = {t-AMl + gli5t-8Ml)} J(i) 



{3t'^ -12tMl + 8M^) g{t) 



4" (t) = {2gl{5Ml-2t) + {l + Ac,m^)it-A Ml) } J(t) 

+ JA!I^(t-4Ml)g{t) 
2F^M^^ "^^^^ 

The functions Df(t) and -62^ (t) are then obtained by subtracting from the 
above expressions the first few terms of their expansion in t: 

Dtit) = 4o + 4i t + 4i t' + Dtit) 
D3it) = dl, + dl,t + D^it) 
B2it) = hl, + B^{t) 

K Integral equations for ttAT scattering 



In section ^ we briefiy described a system of integral equations for the 
partial waves of vrA^ scattering, which is analogous to the Roy equations for 
the TTTT scattering amplitude. The present appendix provides the technical 
details needed to explicitly work out the kernels that occur in these equations. 

Our starting point is the hypothesis that, at low energies, the properties 
of the scattering amplitude are governed by the pole from the Born term and 
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by the unitarity cuts in the s-, t- and w-channels and that only the S- and 
P-waves generate a significant contribution to the latter. In the s-channel 
partial wave decomposition, the terms from these waves read 

\ E + mr^ E — rriff 

= An SStllllk + f^^ljt + . . . 1 , (K.l) 

where z = cos 6 is related to the momentum transfer, t = 2q^{z — 1). The 
variables q and E denote the momentum and the energy of the nucleon in 
the centre of mass system, respectively. The analogous terms from the lowest 
t-channel partial waves read 

^+ ^ 167r/° ^ ^„ ^ 247rm^z/(v/2m^/i -2/1) ^ 

4m^ — t ' 4m^ — t ' 

5+ = + ... , fi- = 6V2 7r/i + ... , (K.2) 

where the dots stand for contributions with J > 2. 

This shows that the relevant contributions from the singularities in the 
s- and M-channels are linear in t, while those from the t-channel are linear 
in V. The above hypothesis is thus equivalent to the assumption that the 
invariant amplitudes and may be described in terms of functions of 
a single variable, which moreover only have a right hand cut:[] 

= A;, + At{s) + At{u) + tAt{s)+tA+{u) + At{t) + A; , 
A- =A-^ + A-{{s)-A-{{u)+tA-:,{s)-tA-:,{u) + uA-^{t) + A~, (K.3) 
B+ = fi+ + Bt{s) - Bt{u) + t B+{s) - t B+{u) + B^ , 
B- = B;^, + B^{s) + B^{u)+tB,{s)+tB,{u) + B^{t) + B; . 

The first terms on the right denote the pseudovector Born contributions, 
while A^ and B^ stand for polynomials in the Mandelstam variables. The 
imaginary parts of the functions A^{s), . . . , B^ (t) are determined by those 
of the 9 partial waves occurring in ( |K.1|) , (|K.2| ): In the notation of eqs. ( P3.3|) , 



^In the present context, the invariant amphtudes are more convenient to work with 
than D^, because the P-waves /{^ give rise to a quadratic t-dependence in D^. 
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(|23.4|), the discontinuities of Af (s) , Af (s) , Bf (s) , Bf (s) are given by 



T .± ^ . , y/j + m^.) (Im/f + 3 Im/3^) O-^ft - ^^ft] 

1 ZTi I ™ 



E + m^r E — rrij^ 

Im^t = 6. 'f + (K.4) 

+ mi^Y[E - rriM) 

ImBf = 4n — — + Air ^ '''' 



E + 771^ E — mj. 



ImBt = Gtt - — 

2 (E + m^)2(E-m, 



while those of A^ (t), B^ (t) are determined by the t-channel partial waves: 
Im/4+ V2m^Im/4--2Im/5- 



ImA+ = 167r , ImAj = -247rm„ , o . 

4m^ — t 4m^ — t 

ImB^ = , Im53- = 6V2 7ilmf^ . (K.5) 

We can now set up the analogon of the Roy equations. In order not to 
burden the discussion with problems of technical nature, we first disregard 
subtractions. The dispersion relations obeyed by the functions associated 
with the s-channel cuts then read 

G{s) = i r d^'^^Gis') ^^^^^ 

J(mjv+Af,)2 S' - S-ie 

while those related to the t-channel singularities are of the form 

^ J4M2 V -t-ie 

The above equations specify the amplitude in the low energy region, in terms 
of the imaginary parts of the 9 relevant partial waves and of the constants 
occurring in the polynomial parts A^.B^ of the representation ( [K.3| ). In 
particular, the real parts of those partial waves may be worked out from the 
s- and t-channel partial wave decompositions of that representation. 

To extract the s-channel partial waves, we need to perform the integrals 



/+1 /•+! 
dzA^{s,t)Pt{z), I3f{s) = j dzB^{s,t)Pe{z) , (K.8) 

where t is to be expressed in terms of s and z, with t = 2q'^{z — 1). The 
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partial waves of interest are then given by 

+^^{-af{s) + {^s + m^)f3t{s)} , 

ftis) = ^^{af{s) + {V~s-m^)PHs)} (K.9) 
^"'""{-4(s) + (Vi + m.)/3o^(s)} , 



ftis) = f^K(«) + (V^-m.)/?^(.)} 

+ {-4(g) + {y/^ + m^) Piis)} . 

In order to extract the i-channel partial waves, we need to expand the 
amplitude in powers of p at fixed t. Concerning the contribution from the 
i-channel cuts, the operation is trivial, because these are linear in u. Those 
arising from the cuts in the s-channel represent a superposition of terms of 
the form 



s' -s s' -ml-Ml- 2mjv i' + \ t 



with s' > {rriN + M^^y. The relevant angular variable is given by 



= - cos^t , p_ = \Jml-\t , q- = ^jMl-\t , 



so that the partial wave expansion amounts to the series 
1 1 ^ 



5^(2£+l)P,(Z)Q,(Z') , (K.IO) 



s'-s '^P-q- 

^, ^ 2s' - 2ml -'^Ml + t 

In the Z— plane, the expansion converges in an ellipse that has the foci at 
Z — ±1 and passes through the point Z — Z'. In the context of our system 
of integral equations, we need the channel partial waves for t > 4M^. For 
real values of Z and t < 4m^ , the variable u is then purely imaginary and the 
denominator s' ~ s is different from zero, in the entire range of integration 
over s'. Hence the partial wave expansion converges and can be inverted 
without problems. 
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In fact, in the case of the amphtude C, we may apply the formula ( [A .21 ) 
to convert the dispersive representation into a series of Legendre polynomials 
of the form ( [A. 21 ) and read off the coefficients. For the amplitude B, however, 
the partial wave decomposition involves the first derivative of the Legendre 
polynomials. The corrolary of ( K.lOj ) that is relevant in this case reads: 

Vp;(Z){Q,_i(Z') -Qmrni • (K.ll) 



1=1 



This then allows us to read off the contributions to t— channel partial waves 
that are generated by the s- and w-channel cuts. 

For the specific partial waves of interest, we may also write down explicit 
projection formulae analogous to (|K.8|) and (|K.9|) : 



57r 



V2 r 



flit) = J7- dZ{l-Z')B-{u,t), (K.12) 



-1 



Wit) 



1 

/ dZZ{p_A-{iy,t) + q_m^ZB-{iy,t)] 

87rg_ y_i ' 



with z/ = Zp_g„/mjv. Note that, for 4M^ < t < 4m^, the integral extends 
over purely imaginary values of u. 

Collecting the various pieces, the system of equations takes the form 
( |23.5| ): The S- and P- waves are expressed as linear superpositions of contri- 
butions from the Born term, from the polynomial and from the imaginary 
parts of these waves. The advantage of working in the t-channel isospin ba- 
sis is that the system does not intertwine the partial waves with even and 
odd isospin. For the t-channel partial waves, the unitarity condition is also 
diagonal in these variables: 

Im/+(t) = {1 - AM'Jty^^ t',{tr ft{t) , (K.13) 
Im/-(t) = {l-AMl/t}-lt\{tyf-{t), z = 4,5 

In the s-channel, however, unitarity does connect the two sets of amplitudes. 
The relation ( 14. 4| ) applies to the partial waves of a given s-channel isospin: 

f} = f^ + 2fr, f} = f+-fr, z = l,2,3. (K.14) 

Expressed in terms of the amplitudes f^, f^, elastic unitarity requires 

Im/+(s) = q\ftis)\' + 2q\fr(s)\', (K.15) 
Im/-(.) = qf+{s)fr{sy + qf7{s)f^{sy + q\fr{s)\'. 
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Finally, we need to specify the polynomial part of the representation 
(|K.3|) , which accounts for the contributions from the higher partial waves. 



This part is intimately connected with the issue of subtractions, because the 
contributions from the subtraction constants amount to a polynomial. Note 
also that the decomposition of the amplitude in eq. ( |K.3|) is not unique. We 
may, for instance, add a constant to the functions A'^ (s) , A2 (s) , (s) , (s) , 
without changing the scattering amplitude. 

While the imaginary parts of the functions Bf(s) fall off sufficiently fast 
at high energies for the integrals in eq. ( [K.6|) to converge, the dispersion 



relations for the functions Af(s) need to be subtracted at least once. In order 
to suppress the contributions of the higher partial waves we oversubtract and 
introduce two subtractions for the amplitudes Af{s) and one for Bf{s). omce 
the contributions of Af{s) and Bf{s) are suppressed by a factor of t, we use 
the minimal number of subtractions for these. Crossing symmetry reduces 
the number of subtraction constants by a factor of two, so that the eight 
subtractions lead to four subtraction constants. We identify these with the 
subthreshold parameters d^Q, diQ, d^^, Bqq, so that we can write the polynomial 
part of the amplitude as 

A+ = <o + t<, b; = o (K.16) 

=^ (doo - &00) ' Bp = , 

Since the amplitude B^ is odd under crossing, we can perform an addi- 
tional subtraction in without introducing a new subtraction constant. 
Subtracting at Sq = rn^ + M^, the dispersion relations in the variable s take 
the form 

Gn{s) = / )n(,f_,) ' (^-1^) 

^ J{mM+M„)2 [S - So) [S - S) 

with n = 2 for Af[s) and n = 1 for (s), -B^(s), Only the one for 

remains unsubtracted, n = 0. The analogous dispersion relations in 
the variable t read 



with n = 2 for A^{t) and n = 1 for A^{t), B^{t). 

This then completes our system of equations. In order to solve it, the 
system must be iterated. Starting with a given input for Im/;^, the rela- 
tions ( |23.5| ) determine the corresponding real parts. The unitarity conditions 
( |K.13D , ( [K.15D then yield a new set of imaginary parts, etc. Note that the 
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outcome is meaningful only at low energies, whereas the dispersion integrals 
extend to infinity. Also, elastic unitarity is valid only at low energies - we 
need experimental information for the behaviour of the imaginary parts of 
the partial waves in the inelastic region, s > {rrij^ + 2M^)^, t > 16M^. 

In the case of vrvr scattering, the mathematics of the Roy equations has 
been explored in detail and it has been shown that two subtractions suffice 
to arrive at a remarkably accurate representation in the elastic region 



For vrA^ scattering, the properties of the corresponding system of equations 
yet need to be explored. Concerning the behaviour of the imaginary parts 
at higher energies, as well as the contributions from higher partial waves, 
the phenomenological information is much better than in the case of vrvr 
scattering. Also, the tttt phase shifts, which enter the t-channel unitarity 
condition, are now known very accurately [0. We are confident that the 
proposed system of equations will indeed provide the required link between 
the physical region and the Cheng-Dashen point and thus allow an accurate 
determination of the cr-term. 
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